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Abstract 

Since Wilson's work on lattice gauge theory in the 1970s, discrete versions of field theo- 
ries have played a vital role in fundamental physics. But there is recent interest in certain 
higher dimensional analogues of gauge theory, such as p-form electromagnetism, including the 
Kalb-Ramond field in string theory, and its nonabelian generalizations. It is desirable to dis- 
cretize such 'higher gauge theories' in a way analogous to lattice gauge theory, but with the 
fundamental geometric structures in the discretization boosted in dimension. As a step toward 
studying discrete versions of more general higher gauge theories, we consider the case of p-form 
electromagnetism. We show that discrete p-form electromagnetism admits a simple algebraic 
description in terms of chain complexes of abelian groups. Moreover, the model allows discrete 
spacetimes with quite general geometry, in contrast to the regular cubical lattices usually as- 
sociated with lattice gauge theory. After constructing a suitable model of discrete spacetime 
for p-form electromagnetism, we quantize the theory using the Euclidean path integral formal- 
ism. The main result is a description of p-form electromagnetism as a 'chain field theory' — 
a theory analogous to topological quantum field theory, but with chain complexes replacing 
manifolds. This, in particular, gives a notion of time evolution from one 'spacelike slice' of 
discrete spacetime to another. 



*A version with better graphics and some related materials can be obtained at http://math.ucr. edu/^derek/pform 
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1 Introduction 



For the relativist, describing electromagnetism as a gauge theory has obvious appeal. The 
language of fiber bundles in which gauge theory is written provides a coherent framework for 
doing field theory in curved spacetimes of arbitrary dimension and with very general global 
topology. The bundle formalism emphasizes the idea that there is no canonical way to compare 
'states at x' with 'states at y' when x and y are distinct points of spacetime — the notion of a 
connection is central. Briefly, gauge theory has much of the sort of flexibility that relativists 
like to demand. 

In contrast, discrete analogs of gauge theory tend to be much more rigid. Field theories are 
often discretizcd using lattice methods [23], in which Euclidean spacetime is approximated by 
a regular hypercubical lattice aZ^, where a is some fixed length called the lattice spacing. The 
lattice way of doing field theory has serious computational advantages, particularly because 
computer models are by nature discrete. But lattice gauge theory also represents a certain 
theoretical compromise. As A. Zee has pointed out (see Chapter VII. 1 of his textbook [29]) 
lattice gauge theory rescues gauge invariance from the "mangling" it suffers in perturbative 
field theory — including the introduction of unphysical 'ghost fields' necessary to make sense 
of the results — but does so at the expense of Lorentz invariance. The discretization shatters 
the symmetry of spacetime. The lattice picks out preferred directions. And the situation is 
worse for those interested in spacetimes more general than afRne Minkowski space: there is 
no such thing as a diffeomorphism-invariant lattice. Of course, lattice gauge theory is often 
considered merely as a computational tool for doing calculations in ordinary continuum gauge 
theory. One expects to recover Lorentz invariance in the limit a — > 0, after Wick rotation. 
If we are really interested in manifestly discrete theories of physics, however, the success of 
lattice gauge theory in its usual form is not very satisfying: we would like to have 'unmangled' 
gauge invariance without doing such drastic damage to spacetime's symmetry. In the context 
of quantum gravity, where spacetime itself is a dynamical variable, the 'spin foam' approach 
avoids picking out preferred directions in spacetime by integrating over possible discretizations; 
in lattice electromagnetism we should at least like to have a theory which allows discretizations 
more general than the hypercubical lattices of lattice gauge theory, with the hope of recovering 
something that looks like diffeomorphism invariance at large scales. 

More generally, we would like such a flexible discretization of p-form electromagnetism, 
the generalization of Maxwell's theory one obtains by replacing the gauge field Ahy a p-form 
[7]. More precisely, the gauge field in electromagnetism is a connection on a C/(l)-bundle, 
which may be thought of in local coordinates as a 1-form 

.A — A-jj^dx^ ^ 

and this may be generalized to a p-ioim: 

A = A^,,„f,^^dx^'' A ■ ■ ■ A dx"" . 

The electromagnetic field, locally the 2-form F = dA, is globally just the curvature of the 
connection A. In direct analogy, the p-form electromagnetic field is described locally by the 
{p + l)-form 

F = dA. 
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The classical equations of motion for p-iorm electromagnetism, what we might call the p-form 
Maxwell equations, look formally identical to the 1-form case [3]: 

dF = 
*d*F = J 

But here, since is a (p + l)-form, -kd-kF is a p-form, so the current J must also be a p-form. 
One sees how the dimensions get boosted: whereas the current in ordinary electromagnetism 
is fundamentally 1-dimensional, corresponding to 0-dimensional point particles tracing out 1- 
dimensional curves, 2-form electromagnetism has a 2-dimensional current. One might expect 
such a theory to describe 1-dimensional charges sweeping out surfaces in spacetime. It is 
thus not surprising that something like this shows up in string theory, where it is called the 
'Kalb-Ramond field' [13, 18]. Integrating the Kalb-Ramond field over the string's worldsheet 
gives a term in the action, analogous to the term in the action of a charged particle given by 
integrating the vector potential A along the particle's worldline. 

A natural question which has been pursued by Baez and others is how to describe p- 
form electromagnetism, and other 'higher gauge theories' [4, 25] with more general gauge 
group, globally. In fact, putting 'p-form gauge theory' on the same global footing as ordinary 
gauge theory requires some higher-dimensional generalization of the notion of fiber bundle 
[4, 5, 6, 19]. In particular, one needs to know what is meant by parallel transporting a p- 
dimensional extended object along a (p + l)-dimensional submanifold, rather than merely 
translating point particles along paths. 

Remarkably, many of the technical difficulties of defining 'higher gauge theory' melt away 
in a discrete setting when the gauge group is abelian, that is, in the case of discrete p-form elec- 
tromagnetism! The need for the bundle formalism in gauge theory arises from the possibility 
that the gauge field might be 'twisted' in a global way. But in the lattice context, changes in 
the gauge field happen not continuously but in discrete steps. It therefore becomes impossible 
to decide if the field is 'twisted' around some noncontractible loop. Moreover, the abelian 
nature of electromagnetism allows higher-dimensional generalization to be carried out in a 
straightforward way. Rather than increasingly rich algebraic structures necessary for boosting 
the dimension in nonabelian gauge theory, higher dimensional electromagnetism can always 
be described in terms of chain complexes of abelian groups. In fact, p-form electromagnetism 
on discrete spacetime turns out to be no more difficult than ordinary electromagnetism. 

It is the hope of the author that this detailed and pedagogical study of discrete p-form 
electromagnetism will serve as a starting point for further study in both lattice gauge theory 
and higher-dimensional analogs of gauge theories, both fields of active research, as well as 
theories at their interface. A subsidiary motivation is to smooth the reader's way toward an 
understanding of discrete models of spacetime which arise in the study of quantum gravity, 
especially 'spin foams'. It is also hoped that the expository flavor of this article will make it 
accessible to a wider variety of readers, including physicists wishing to see an elegant treatment 
of lattice gauge theory and mathematicians wishing to see applications of cohomology and 
category theory in physics. 

The plan of the paper is as follows. Section 2 serves as an invitation to much of the 
material that follows it, with the intent of getting as quickly as possible to the point of doing 
path integral calculations in discrete quantum electromagnetism. In this section we make 
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an initial attempt at quantizing discrete vacuum electromagnetism with gauge group M. As 
typical in lattice gauge theory, we work in this section and throughout the paper in Euclidean 
(i.e. 'Wick-rotated') spacetime. Readers already well acquainted with path integrals in lattice 
gauge theory may wish to skim through this to fix notation and then move on to Section 3. 

In Section 3 we begin a more careful approach, considering first what mathematical model 
of discrete spacetime is appropriate for the description of p-form electromagnetism. 

Section 4 draws a comparison to the continuum theory and uses this to construct an 
appropriate form of the action in discrete p-form electromagnetism. In a subsection, it is 
shown that, just as ordinary gauge theories are particularly simple in two dimensions, p-form 
electromagnetism simplifies greatly in p + 1 dimensions. 

Topological issues are considered in Section 6, where the cohomology of a spacetime is 
introduced as a criterion for convergence of the path integral when the gauge group is M. The 
p-form analog of the Bohm-Aharonov eflFect is also discussed. In the final subsection of Section 
6, the classical theory is examined more closely, using the topological machinery developed 
earlier in the section. Using Hodge's theorem, it is shown that, quite surprisingly, under 
certain conditions p-form electromagnetism in Riemxmnian signature is purely topological as 
a classical theory. More precisely, in the absence of charged matter the space of classical 
solutions depends only on the topology of spacetime. 

In Section 7 we examine the relationship between the groups M and U{1) as gauge groups 
for electromagnetism, and use the formalism developed in previous sections for gauge group 
M to set up path integrals in the more usual U{1) case. These these turn out to involve not 
Gaussian integrals in the ordinary sense, but the analog of a Gaussian on a torus, which is 
nothing but a theta function. 

Sections 8 and 9 deal with some special cases of the theory. First, we consider 0-form 
electromagnetism, which is really just scalar field theory, and find that understanding this 
requires a further insight in topology, namely, augmentation of complexes. We then return to 
the case of p-form electromagnetism in p + 1 dimensions, a case which is in fact the source 
of many examples in the paper. We suggest that in p + 1 dimensions, the theory is 'almost' 
a topological quantum field theory — what we call a 'volumetric field theory' — where the 
volume of spacetime is the only nontopological degree of freedom. 

The final section of the paper. Section 10, presents the main result of this work, namely a 
description of lattice p-form electromagnetism that parallels the category-theoretic description 
of topological field theory. The importance of this result is that it gives a rigorous description 
of time evolution in lattice p-form electromagnetism. Namely, it allows an n-dimensional 
discrete spacetime to be split into (n — l)-dimensional slices representing space, and provides 
time evolution operators between the Hilbert spaces of states on these. Though some attempt 
has been made to retain the expository tone of the previous sections, this final section is 
more technically demanding than the rest of the paper. In particular it will probably not be 
comprehensible without a working knowledge of categories and functors. A highly readable 
introduction to categories and functors which also emphasizes applications to physics is the 
book by Geroch [11]. The more technical aspects of category theory needed in this section are 
all covered in Mac Lane's textbook [20]. 

As this paper is intended to be expository and as self-contained as possible, an Appendix 
includes a review of the Gaussian integrals needed to do calculations in discrete p-form elec- 
tromagnetism. 
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2 Survey: From n-Graphs to Path Integrals 



A discrete model of spacctime might begin with a simple discrete set of points, the "events." 
Soon however, if we want to write down theories that look like their desired continuum limits, 
we are forced to equip our model with additional data. In the context of gauge theory, for 
example, we would like to be able to define an analog of the connection and take holonomies 
along paths, and this requires that we specify ways of getting from one event to another. This 
has the effect of turning our model into a graph, with events as vertices and paths between 
events as finite strings of e^e 




Likewise, the desire to give meaning to notions such as area and curvatu?5^akes it natural 
to fill in empty spaces bounded by edges with 2-dimensional faces, or plaquettes. Thus our 
model becomes a '2-graph': 




In a similar way, we are led to adjoin 3-dimensional cells bounded by plaquettes, and so on, 
all the way up to the highest dimension of the spacetime. In general, an n-dimensional discrete 
spacetime should be constructed from n+1 distinct classes of geometric objects: vertices, edges 
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between vertices, plaquettes bounded by edges, 3-cells bounded by plaquettes, and n-cells 
bounded by (n — l)-cells. Spacetime thus becomes an 'n-graph'. 

The problem is, nobody knows what an n-graph is! Of course, we have some soHd ideas of 
what an n-graph should be like, but there seems to be no accepted intrinsic definition that is 
flexible enough to encompass every would-be example. We take up this issue in Section 3. 

For now, let us rely on our intuitive notion of a 2-graph as a combinatorial object consisting 
of finite sets of vertices, edges, and plaquettes, and begin with a naive attempt at discretizing 
quantum elcctromagnctism in 2 dimensions. Let V be the set of vertices, E the set of edges, 
and P the set of plaquettes or faces. An edge is just a line segment or curve connecting 
one vertex to another, while a plaquette is a 2-dimensional surface whose boundary consists of 
one or more edges. A simple example, with 2 vertices, 3 edges, and 2 plaquettes, is depicted 




In practice, it is helpful to arbitrarily assign an orientation to each cell. It is also useful 
to give labels to cells, just as it is sometimes convenient to use coordinates in ordinary gauge 
theory. We shall label the elements of V, E, and P using lower case letters v, e, and p, with 
subscripts to distinguish between them, as shown below: 



V2 




Vl 



The usual gauge group for electromagnetism is U{1), but there is another obvious choice, 
namely M. We will follow the latter alternative in the first part of the paper, both because a 
noncompact gauge group leads to some instructive problems, and because the construction of 
the U{1) path-integral in Section 7 is made easier having established the formalism with gauge 
group M. In lattice gauge theory, the gauge field or connection A assigns to each edge in 
the lattice an element of the gauge group: 

A: E^R. 

The space of connections is therefore just 

A = R^. 

Now suppose we move a test particle from one vertex to another along some path, i.e. along 
some concatenation 7 = e^^e|^ . . . e^" of oriented edges, where = ±1 with e^^ denoting the 
edge ei = ef^ with its orientation reversed. We define the holonomy of A along the path 7 
to be the sum 

H{A, 7) = eiA{ei) + £2^(62) + • • • + £n^(en) 
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Note that we use addition and subtraction because R is an additive group. Most gauge 
theories, including the ?7(l)-version of electromagnetism, use multipUcative groups, so the 
holonomy becomes H{A,j) = ^(61)^^^(62)^^ • • • A{enY" instead. 

In gauge theory, one way to tell if a connection is curved is to send a test particle around 
some small loop and keep track of what happens to the connection as it goes around. In the 
context of electromagnetism, this means comparing the phase of the particle before and 
after its trip around the loop. In the present context, our gauge group is not the group of 
phases, but the idea is similar. Accordingly, wc define the curvature or field strength F{p) 
on a plaquctte p G P to be the holonomy of the connection A around the loop that forms its 
boundary dp. For example, in the lattice 



the curvature on plaquette pi is A2—A1, while the curvature onp2 is As—A2, where Ai := A(ei). 
Note that since the gauge group is abelian, the holonomy around the loop does not depend on 
what order we traverse the edges in! This fact will be of great use to us. 

To quantize our theory, we will use the Euclidean path integral approach, and to do this 
we must specify the form of the action for the gauge field. Roughly, when there are no 
sources present, the action should be a positive quadratic function of the curvature, so that 
the connections with zero curvature have the least action, hence are most probable. In Section 
4 we discuss in more detail what the action for a given connection should look like. In the case 
of lattice electromagnetism in two dimensions, we find that the best choice for the action is 



where Fj = F{pi) is the curvature induced by the connection on the ith plaquette, is the 
area of the ith plaquette, and is the square of the charge on the electron (or if you wish, 
the fine structure constant a = e'^/hc, since we are using units with h = c = 1). 
An observable O in our theory is a real-valued function of the gauge field: 



In certain cases, it might be difficult to define exactly what this means. For example, if 
our lattice is infinite, then we must deal with the difficulties of defining measures on infinite 





O: A^R 



so we can try to calculate its expected value, using our action S{A), as 
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dimensional Euclidean space. However, even in the finite dimensional case, where VA = dPA 
is just Lebesgue measure on ^ = M-^, when we try to calculate (O) by the above formula we 
usually run into serious problems! To see this, let us do a sample calculation of the integral in 
the denominator, which is just the partition function 



against which all other expected values are normalized. 

For convenience, let us take the little two-plaquette complex above as our spacetime, with 
plaquette areas Vi and V2. We then write the action as a function of the connection using 
matrix notation: 

2e2 V ^1 ^^^2 



2e^ViV2 
1 

2e'^ViV2 



[Ai A2 A3] 
{A,QA) 



V2 -V2 
-V2 V1 + V2 -Vi 
-Vi Vi 









A2 




A3 



where A in the last line denotes the vector with components (^1,^2,^3) and Q is the 3x3 
matrix in the line above. If we naively apply one of the integal formulas derived in the 
Appendix, we obtain 



-{A,QA),E^ 



d^A 



\2TTe'^ViV2f 
det(Q) 



However, when we actually calculate the determinant we find det((5) = 0, thereby obtaining 
the most famous wrong answer in quantum field theory: Z = oo\ This is obviously bad. 

Let us examine why the integral diverges. Notice that if we transform the connection as 
follows: 

(^1,^2,^3) ^ iAi + d(l),A2 + d(t),A3 + d(f)), 



where d(p is any real number, the action 



SiA) 



1 

2i2 



(^1 



Vi 



+ 



{A2 



V2 



is unaffected. This symmetry of the action is a degree of 'gauge freedom'. From this per- 
spective, the divergence of the partition function should hardly be surprising — the integrand 
exp(— 5) is equal to 1 along an entire line through the origin of R^! More pictorially, one cross 
section of exp(— 5) looks like this: 
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where the direction of gauge freedom is obvious. 

The reason we could add {d<p, dcp, d(j)) to the connection A = (Ai, ^2, ^3) without affecting 
the action is that the connection {d(p, d(p, d(p) is flat: its holonomy around each plaquette is 
zero. A Uttle thought reveals a simple way of constructing such fiat connections. First, for 
each directed edge e in the lattice, let s{e),t{e) G V denote the source and target (or 'tail' 
and 'tip') of e, respectively: 

s( e) e tie ) 

Then, if we let 

^: F^M 

be any real-valued function of the vertices of spacetime, and define the differential of to 

be 

d(t): 

given by 

d(t>{e) = me)) - Hs{e)) 

the connection dcj) is flat! To see this, note that computing the holonomy of d(j) around any 
loop then amounts to just adding and subtracting the same value at each vertex along the 
way. For example taking the holonomy of d(f) along the boundary 7 of this plaquette: 

^3 

62/ \ei _i 

7 = 6162 63 

Vl 63 V2 

gives 

H{d4),-f) = #(ei) - d(^(e2) + #(63) 

= {cp3 - h) - ih - ^1) + ih -cpi)=0 

where (pi := 4>{vi) for i = 1,2,3. In fact, in lattice gauge theory wc define a gauge trans- 
formation to be a map from the set of vertices to the gauge group. The space of gauge 
transformations is thus 

g = M^. 

It is convenient to think of ^ as a group which acts on the space A of connections by 

{(p, A) t-^ A + d(j). 
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A common way of eliminating divergences caused by gauge freedom is gauge fixing. Fixing 
a gauge means choosing some method to pick one representative of each gauge-equivalence class 
[A] and doing path integrals over just these. There are problems with this approach, however. 
In a general gauge theory, it might not even be possible to fix a smooth, global gauge over 
which to integrate. But even when we can do this, the arbitrary choice involved in fixing a 
gauge is undesirable, philosophically. Gauge fixing amounts to pretending a quotient space 
of A, the space of connections modulo gauge transformations, is a subspace. A better 
approach is to use the quotient space directly. Namely, modding out by the action of Q on A 
gives the quotient space A/G consisting of gauge-equivalence classes of connections, and we 
do path integrals like 

JA/g 

When gauge fixing works, path integrals over A/G give the same results as gauge fixing. But 
integrating over A/G is more general and involves no arbitrary choices. 

Since gauge equivalent connections are regarded as physically equivalent, the quotient space 
.4/^ of connections mod gauge transformations is sometimes called the physical configura- 
tion space for vacuum clectromagnetism. A physical observable is then any real-valued 
function on the physical configuration space: 

O-.A/G^'R 

or equivalently, any gauge-invariant function on the space A of connections. 

In the case of noncompact gauge group even factoring out all of the gauge freedom may 
be insufficient to regularize our path integrals. In particular, there are certain topological 
conditions our spacetime must meet for this programme to give convergent path integrals. As 
a consequence, in the case where the gauge group is M, we must sometimes take more drastic 
measures in order to extract meaningful results, and this leads to some interesting differences 
between the cases G = M and G = U{1) as the gauge group for clectromagnetism. These 
differences are related to the famous 'Bohm-Aharonov effect', which has higher dimensional 
generalizations in the case of p-form clectromagnetism. We discuss this issue in Section 6. 

3 Spacetime Lattice as an n-Complex 

Just as the concept of a manifold is fundamental to an elegant treatment of general relativity, 
lattice gauge theory must give due attention to the spacetime lattice itself. In particular, we 
should say with some precision what we mean when we describe discrete spacetime as an 'n- 
graph.' We make no attempt in this section at giving the most general definition of n-graph 
suitable for describing arbitrary discrete gauge theories. Rather, we show that for abelian 
gauge theory, we really only need a weaker notion — what we shall call an 'n-complex'. 

3.1 The n- Graph Problem 

Formulating a precise definition of n-graph is not an easy task. Ideally, we should be able to 
describe n-graphs in a purely combinatorial way: An n-graph should consist of sets of various 
kinds of cells, together with maps telling how the cells are linked together. An ordinary 
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(directed) graph, for example, is specified by a set V of vertices and a set E of edges, together 
with source and target maps 



telhng at which vertex each edge begins and ends. 

If one tries copying this definition to define a notion of n-graph, one quickly realizes the 
immensity of the combinatorial problem at hand! Even if we only wish to add 2-dimensional 
faces, complete generality requires an infinite number of new sets and set maps; whereas an 
edge has only two ends, a 2d face can have arbitrarily many sides: 



What's more, the specified maps must satisfy certain equations. For example, in the pentagonal 
face 



we must impose equations saying that the 0th and 4th edges have the same source vertex, that 
the target of the 2nd edge is the source of the 3rd, and so on (three more equations). These 
equations are called incidence relations. 

And it gets worse. A 3-cell may have arbitrarily many faces, each of which may have any 
number of sides . . . 

There are various attempts at formalizing the notion of n-graph in a fully satisfactory way, 
most relying heavily on category theory, or rather n-category theory, which itself is still under 
development. At the present state of the art, however, we are forced to make compromises of 
one sort or another. Each of these compromises carries with it its own weaknesses. 

We could, for example, restrict the allowable shapes of cells. Allowing only simplices as 
cells, for example, lets us define 'simplicial n-graphs' as simplicial sets [21], where we simply 
remove any cells of dimension greater than n. Similarly, we could define 'cubical n-graphs' 
using cubical sets [12]. If we are only interested in lattice gauge theory as a computational 
tool, these should be sufficient, and indeed the study of lattice gauge theory began with regular 
cubical lattices. However, if we really believe spacctimc is discrete at some scale, then there 
seems to be no a priori physical justification for imposing such strong conditions on the shapes 
of cells! 

For this reason, piecewise-linear CW-complexes (or PLCW-complexes, for short) [16, 17] 
have been used as discrete models of spacetime in loop quantum gravity [1 , 2] . These allow 



s,t: E^V 




• • • 
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a wider variety of cell shapes, and are quite handy if, for example, we want to chop up a 
manifold with boundary into polyhedra: 




Unlike simplicial sets, however, we do not know how to define PLCW complexes intrinsically 
— they rely on an ambient space for their construction. This too is undesirable for physics, 
analogous to defining the spacetime manifold in relativity as a submanifold of some higher 
dimensional Euclidean space. In addition, the polyhedral cells in PLCW complexes would not 
allow, for example, plaquettes bounded by only one or two edges: 




or 3-cells bounded by fewer than four plaquettes, like these: 




and so on. The problem with all of these is that they involve curved sides; cells in PLCW 
complexes must be constructible in Euclidean space using only polyhedra with straight sides. 
So in fact, even our simple two-plaquette example from Section 2 is not a PLCW complex. 
Worse, a typical graph is not a 1-dimensional PLCW complex, so PLCW complexes are hardly 
a candidate for general n-graphs. 
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3.2 n-Complexes 

Luckily, it turns out that in gauge theory, and even 'higher gauge theory', we can avoid 
many of the technical issues with n-graphs when the gauge group is abelian — that is, in the 
case of p-form electromagnetism! To see why, note that taking the holonomy of a connection 
along some path through a graph, we multiply group elements of the edges composing the 
path in order. In a nonabclian gauge theory, this rcqTiircs we keep close track of how the edges 
are interconnected. This is manageable for 1-dimcnsional holonomies, but in "p-ioxva. gauge 
theory', where holonomy involves multiplying group elements labelling p-cells, it is difficult 
to know what order to use. In general, the holonomy depends intimately on the incidence 
relations of the n-graph — that is, on the details of how cells are linked up — and keeping 
track of this requires a fairly sophisticated notion of n-graph. If the gauge group is abelian, 
however, then the order in which we apply the group operation is unimportant, and all that 
matters is the orientation of each cell! 

The key to defining a rudimentary n-graph suitable to our purposes is to follow the topolog- 
ical maxim, 'The boundary of a boundary is zero '. This principle, so phrased by John Archibald 
Wheeler [8, 22], has limitless applications to physics, from Kirchoff's laws in electrical circuits 
to the Bianchi identities in general relativity. It is also the topological foundation of Maxwell's 
equations themselves, so it is perhaps not so surprising that it plays a foundational role in a 
precise formulation of lattice electromagnetism. 

For us, spacetime will be a cell complex M, which we define as follows. First of all, M 
consists of a list of sets 

Xo,Xi, X2 ■ ■ ■ , Xk, . . . 

where we call Xk = Xk{M) the set of fc-cells in M . Now intuitively, the boundary of a /c-cell 
should be a "sum" of {k — l)-cells. To formalize such sums, for each k we let = Ck{M) 
be the free abelian group on the set X^. In other words, Ck just consists of all formal linear 
combinations of /^-cells, with integer coefficients. We call the elements of the fe-chains in 
M. We then hypothesize boundary maps dk- Ck ^ Ck-i and require that these be linear 
over Z. If we let (?o : Co — > be the unique map from Co to the trivial group, then the boundary 
maps fit together like this: 

(7q (J^ ^^2 83 _ _ _ dk Qfc+1 

We will say that the cell complex is n-dimensional, or that it is an n-complex, provided 
Xfc = for all A; > n.^ When the cell complex is n-dimensional, we might as well truncate the 
chain and write simply 

With this construction. Wheeler's principle can be stated in a concise way: 

dk-idkc = 0, 

for every c & Ck- Better yet, when it is clear which boundary maps we are talking about, we 
often drop the subscripts and just write d for all of them. If we also drop the explicit reference 



^Note that X„ can be empty. Thus an n-complex is trivially also an (n -|- l)-complex. 
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to a /c-chain c, we can write simply 

dd = o 

or even 

d^ = 0. 

So we conclude our definition of a cell complex by demanding that the boundary of the bound- 
ary of any fc-chain is zero, in the sense just described. 

In homological algebra, such a sequence of maps d: Ck ^ Ck-i satisfying dd = is called 
a chain complex. This lets us summarize the definition in a concise way: 
Definition 1 An n-dimensional cell complex, or n-complex M is a chain complex 

of free abelian groups Ck{M), each equipped with a preferred basis Xk{M). 

In a certain sense, this definition is perhaps too inclusive. We could take, for instance, the 
sets Xk of cells to be any sets whatsoever and let each boundary map d be the zero map. One 
can also concoct other sorts of perverse examples, all falling squarely within the bounds of the 
definition, but with seemingly no connection between the boundary maps and the geometric 
notion of boundary. 

The point is that, although the result might be strange or impractical, one actually can 
write down a theory of electromagnctism on such n-complcxcs. Our attitiidc will be simply 
that it is up to the user of the definition to decide what additional conditions to impose to 
guarantee more 'reasonable' n-complexes. But to see that our definition of an n-complex at 
least does include the sort of discrete spacetimes we really are interested in, and to show how 
the definition might be used in practice, an example or two might be in order. 

Example: In Section 2, we had the following spacetime: 

V2 




Vl 



But now, instead of simply talking about the sets Xq = V, Xi = E, and X2 = P of vertices, 
edges and plaquettes, we make them into free abelian groups 

Co = Ci ^ iF, and C2 ^ if, 

with ordered bases F, P, respectively. The boundary map d: Ci Co we define by 
assigning to each edge e G Ci its target minus its source: 

d{e) = t{e) - s{e). 
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In particular, since all three edges in this example go from vi to V2, we have 



d{ei) 
d{e2) 
dies) 



V2 — Vl 
V2 — Vl 
V2 — Vl. 



Similarly, we define the boundary of a plaquette p G C2 by taking the sum of the edges 
around its physical boundary, but with opposite sign if an edge points opposite the plaquette 
orientation. Thus, in our example, 

d{pi) =61-62 

9{P2) = 62 - 63. 

It is often convenient to write the boundary maps in matrix form. Relative to the ordered 
bases {fi, ^2} of Cq, {ei, 62, 63} of Ci, and {pi,P2} of C2, they are represented in our example 
by the matrices: 

1 



di 



-1 -1 
1 1 



-1 





1 

-1 



Obviously 8182 is just the 2x2 zero matrix, as required. 
Example: As a second example, take the solid tetrahedron: 




which has four vertices, six edges, four plaquettes, and one 3-cell. Here we have left out 
the labels on plaquettes, to avoid cluttering the diagram, but for convenience we take the 
zth plaquette Pi to be the triangular face whose vertices include all vertices but Vi. On each 
plaquette, let us take the orientation such that its three edges are traversed in numerically 
increasing order (or a cyclic permutation). To specify the orientation of the 3-cell c, we call 
each of its faces positively oriented if its orientation is counterclockwise when viewed from 
inside c. The resulting chain complex then looks like this: 







Z' 



6 .92_ ^4 ^ 



Z 
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with the boundary maps represented by 



the matrices 



= [0 
















-1 


-1 





-1 





1 


-1 








-1 





1 


1 





1 





-1 





-1 


-1 











1 





1 








1 





-1 





-1 








1 


1 








-1 











1 


1 


1 



Again, we can easily verify that dodi, 8182, and ^293 are all zero, so we really do get a chain 
complex. 

Notice that in passing from the physical model of a spacetime lattice to its corresponding 
n-complex, one is actually throwing away quite a lot of information! Given a chain complex, 
it may be impossible to reconstruct every detail of how cells in the original model might have 
been hooked together. Applying the boundary map to a fc-ccll only tells us 'how many times' 
each (k — l)-cell is included in its boundary, with possible cancellation. For example, in the 
cylinder with one plaquette p: 




we have (9p = ei — 62 + 63 + 62 = 61 + 63, leaving no trace of the fact that 62 is actually included 
twice in the physical boundary of p, with opposite orientation. 

We must emphasize that we are free to discard all of these incidence data only because 
electromagnetism is abelian. As we shall soon see, almost everything we need for electromag- 
netism is contained in the chain complex of spacetime. We can even think of discrete spacetime 
as being its chain complex, a perspective which we explore more fully in Section 10. 

One immediate generalization of the notion of cell complex is sometimes useful. Although 
integer linear combinations of cells are perhaps the most sensible from a purely combinatorial 
perspective, we may instead take real coefficients. The then become real vector spaces, 
rather than mere abelian groups, allowing arbitrary real-linear combinations of cells in the 
spacetime lattice. This is useful because it allows us to exploit the nice properties of the field 
M and vector spaces over it in certain computations. More generally, we could take coefficients 
in any ring R. Formal i?-linear combinations of elements of give the free J?-module on 
X^. When we need to be specific about what ring we are using, we write 

CkiM;R) 

for the i?-module of /c-chains in M. 
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4 Discrete p-Connections and the Action 



4.1 Cochains as Forms 

The local theory of continuum electromagnetism is most beautifully written in the language 
of differential forms [3]. In p-form electromagnetism, the analog of the connection A is locally 
a p-form, making the electromagnetic field F=dA a (p+ l)-form. In lattice jo-form electromag- 
netism, we thus define a p-connection to be a map assigning to each p-cell an element of the 
gauge group: 

A: Xp{M) -^R. 

But in the framework introduced in the previous section, we can now speak not only of p-cells 
but also formal linear combinations of them. The real beauty of this viewpoint is we can 
extend A linearly to obtain a linear functional 

A: Cp{M;R) ^ M. 

That is, we can consider a p-connection to be an element of the dual space 

C^'(M;M) := Cp(M;M)*, 

which we dub the space of p-cochains. Likewise, the curvature is a linear functional on the 
space of {p + l)-cells, or an element of CP~^^{M; M). In general, we define 

C^{M-R) ■-CkiM;R)* 

to be the i?-module dual to Ck{M; R) — i.e. the module of all maps Ck{M; R) ^ R that are 
linear over R. When the ring R is clear from the context (typically M or Z), we may simply 
write C'' for C^{M;R). 

To make our discrete theory look as much as possible like its continuous counterpart, we 
need to understand what takes the place of the differential. This construction is standard in 
algebraic topology [14, 15]. Since the boundary maps form a chain complex, 

n do ^ di ^ 3-2 ^83 dn ^ 
U < Oo < Oi < O2 < • • • < C/„ 

we get a chain complex of the dual spaces in the opposite direction, the cochain complex, 
by taking the adjoint of each boundary map: 

^0 dp ^ di ^ ^2 d,2 ^ _ _ _ dn-l dn^ Q 

That is, given u G C", we define diOJ G C"+^ by diUj{x) = u{di-^-ix). Following the language of 
differential forms, we call the maps di differentials. Again, dropping the subscripts on both 
the boundary maps and the differentials usually causes no confusion and we can write more 
cleanly 

df{x) = f{dx). 

This definition justifies the name cochain complex: it is easy to see that dd = implies dd = 0. 

Exploiting the strong parallel between our cochain complex and the DeRham complex of 
differential forms, we will slip between the two vocabularies at will, calling cochains 'forms' 
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whenever convenient. Also borrowing from differential geometry, we will say that a /c-cochain 
uj is closed if do; = or exact if a; = dr/ for some {k — l)-cochain 77. The equation = 
then has its usual DeRham translation, "all exact forms arc closed." 

Differential forms live to be integrated, so we should understand what it means to integrate 
a discrete differential form. There is a beautiful way of thinking about integrating a cochain 
that will clarify the correspondence between the equation dj{x) = f{dx) and the continuum 
theory. The idea is to think of 'integrating' a A;-cochain u over a 'region' (a A;-chain) R as 
evaluating u{R). If we introduce the notation 

JR 

then the definition of the differential d as the adjoint of the boundary map d is just 

duj = u 
JR JdR 

— a discrete version of Stokes' theorem on differential forms! 



4.2 Action 

In ordinary continuum electromagnetism, the Lagrangian for the free electromagnetic field 
F = dA is given by^ 

This equation holds for free p-form electromagnetism as well, where F is the p-form field 
F = dA. To write down a lattice version of this theory, we thus seem to need discrete analogs 
of the wedge product A and the Hodge star operator The Hodge star in differential geometry 
turns g-forms into (n — g)-forms, so we seem to need a map: 

for each gr = 0, 1, . . . , n. Likewise, for an analog of the wedge product of differential forms we 
would want 

A: C« X ^ C«+^ 

Defining these maps seems to be the main obstacle to completely formulating a discrete 
analog of the theory of differential forms. In certain special cases, it is known what these 
operations should look like. For example, for simplicial complexes, the best approximation to 
the wedge product of differential forms is the 'cup product' U of cochains [14]. For the Hodge 
star, a cubical complex is more convenient, since Hodge duaUty is related to Poincare duality, 
and the Poincare dual of a cubical lattice is also cubical.^ 



^In many books, one sees the free eleetromagnetic Lagrangian written as A or in index notation as 
jFfj^yF^^" , without the coupling constant (1/e^). This means either that the authors have hidden the coupling 
constant elsewhere, usually by defining the covariant derivative to be = — ieA^, or that they are interested 
only in the classical theory, where constant multiples of the action have no effect. After quantization, the magnitude 
of e scales the probability of large vacuum fluctuations. 

^For more on discrete differential geometry, we refer the reader to the papers by Forgy and Schreiber [9] and 
Forman [10], and to additional references which they cite. 
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But we would rather not impose any restrictions on the shapes of cells in our spacetime 
model. Forcing the lattice to be too regular goes against relativity's lesson that spacetime, 
at least on the macroscopic level, does not have preferred directions the way a crystal lattice 
does. 

Luckily, if we resist the temptation to dissect the formula for the action 



we can instead discretize the whole thing in one fell swoop! The key observation is that 



defines an inner product, the Hodge inner product, on differential (p+ l)-forms. So in the 
lattice theory, if we endow C^+^ with an inner product 



then we can define the action corresponding to the field F to be 



What this inner product should look like for a given lattice depends on the structure and 
the dimension of the lattice. We shall not dig too deeply into this issue, but just mention a 
few 'reasonable' requirements. For this it is useful to write the inner product as: 



where the implicit sum is over all pairs of {p + l)-cells Xi,Xj. One sensible condition that 
might be imposed on this inner product is a 'locality condition.' That is, we do not expect 
the action to have correlations between distant points in spacetime, so we might require that 
h^^ vanish whenever Xi,yi are 'nonadjacent' (p + l)-cells. Also, when the details of lattice 
geometry are carefully accounted for, we might expect the correlation between two cells to 
depend on the physical angle between those cells. In particular, on a regular cubical lattice, 
we might well expect k^^ to be a diagonal matrix. Note that neither the locality condition nor 
any geometric data used to build the inner product can be derived from the chain complex 
itself. Both constitute additional structure in the spactime model. 

For the sake of simplicity, many of our examples in this paper involve electromagnetism in 
two dimensions. For this reason, wc next turn briefly away from our more general development 
of the theory to consider this special case. More generally, we consider the case of p-iorm 
electromagnetism in p + 1 dimensions and find that the inner product on p-cochains has a 
natural and simple form. 

4.3 The Case of p + 1 Dimensions 

Ordinary gauge theories in two dimensions are famously well-behaved [28]. One of the 
reasons electromagnetism in two dimension is especially nice is that the electromagnetic field. 




{F, G) := F A-kG 



{F,G)=h'^F{xi)G{xj) 
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being a 2-form, can be written simply as a scalar field times the volume form (or perhaps more 
properly, the 'area form'): 



This allows us to choose the best possible inner product on 1-cochains — an exact analog 
of the Hodge inner product — without having to define a discrete Hodge star operator in 
detail. In particular, the only Hodge duals we really need to consider are the trivial ones: 
★vol = 1 and = vol. More generally, and for the same reason, we can do discrete p-form 
electromagentism and avoid any unpleasant entanglements with the Hodge star as long as we 
work in (p + l)-dimensional spacetime. 

To do path integrals, there are several reasonable choices we could make for the action. 
We would like, however, to arrive at a discrete theory that looks as much as possible like the 
more familiar theory of electromagnetism on a manifold. With this motivation, we make an 
unapologetic appeal to the known continuum theory to decide which action to use. 

Since the usual Lagrangian for vacuum electromagnetism is L = -^F A ★F, we find, in the 
two-dimensional case: 



Now, if we want to discretize this theory, say for the purpose of modeling on a computer, 
we clearly want a fine enough lattice that the value of the scalar curvature field / in the above 
equations is essentially constant over any plaquette. If we make this assumption, then the 
action associated with a plaquette pi whose curvature is F should be 



where Vi is the area of the plaquette pi. Similarly, the most obvious way of discretizing the 
curvature is to assign to each plaquette the area integral of its continuum curvature. 



F = f vol. 



L=^{f yo\) A ^{f vol) 





Pi 





Pi 




Pi 



f-Vi. 
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Combining these last two results we find a candidate for the action which now uses only 
variables that are available to us directly in the lattice model! 



To get the total action, we just sum over plaquettes. Recalling from Section 4.2 that the action 
is given by 

the corresponding inner product on 2-cochains for electromagnetism in 2 dimensions is thus 

Fid 



{F,G)= ^ 



PieX2 

The preceding derivation of the action did not rely explicitly on spacetime being two 
dimensional but only on the fact that the curvature is a scalar multiple of the volume form. 
More generally, this happens for p-form electromagnetism whenever spacetime is (p + 1)- 
dimensional. Indeed, in this case we can repeat the above calculations to obtain the obvious 
generalization of the action. The inner product on (p+l)-cochains for p-ioim electromagnetism 
in p + 1 dimensions is thus 

(^'^)= E ^- 

Here Vi := vol{xi), where we now use vol to denote the discrete volume form — a {p + 1)- 

co chain 

vol: Cp+i R 
that assigns to each {p + l)-cell its volume. 

5 Discrete p-form Maxwell Equations 

As with any field theory, we may obtain the classical equations of motion — in this case the 
discrete analog of Maxwell's equations, or rather their p-form generalization — by extremizing 
the action. Formally, the action in our theory is the map 

given by 

SiA) = ^{dA,dA) 

where {—,—) is the inner product on {p + l)-cochains introduced in the previous section. The 
extrema of this function are the ^-connections A such that the differential 



6S: TaCP T, 



s{Ay 
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is the zero map. We get: 



= SS = ^5{dA, dA) 
= i.{6dA,dA) 
= ^(ddA,dA) 



\{SA,dUA) 



e2 



where equahty in the second line follows from the product rule together with the symmetry of 
the inner product, and d"^ is the Hilbert space adjoint 

rft : CP+^ CP 

of the linear operator d. At an extremum, 5S = must hold for all values of the variation 5A. 
So, by our above calculation, 

dUA = 

or 

d^F = 

which is our discrete analog of the vacuum Maxwell equations^ 

-kdi^F = 0. 

It is worth noting the relationship between the continuum and lattice versions of the 
Maxwell equations. In particular, this means understanding the relationship between the 
two kinds of duality showing up in these equations. In Riemannian signature, the Hodge op- 
erator satisfies -k'^uj = {—1)p^'"~p^uj when acting on any p-form cj. So given a p-form A and 
{p + l)-form G, we have 

{dA, G) = j dAM^G 

= (-if+^y AAdi^G 

= {A, *d*G) 
where in the second step we did an integration by parts. Thus 

In fact, the vacuum ;j-form Maxwell equations: 

d^F = 



^This equation corresponds to the two nontrivial vacuum Maxwell equations, which say V-E = and VxB 
The other two equations come from the tautology dF = 0, just as with differential forms. 
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can be simplified further in the present context. If F = dA satisfies this equation, then we 
have: 

= {d^F, A) = {F, dA) = {F, F). 

In Lorentzian physics, we draw no strong conclusions from this calculation. But in the Rie- 

mannian case, the inner product ( , ) is positive definite, and hence the classical equations 
of motion for free discrete p-form electromagnetism reduce to: 

F = 0\ 

This is shocking if we are accustomed to electromagnetism on Minkowski spacetime. It means, 
in particular, that our theory has no analog of electromagnetic waves propagating through a 
vacuum! Moreover, the same proof works not only in the discrete case but also, for example, 
for p-form electromagnetism on a compact Riemannian manifold. The two assumptions leading 
to the conclusion that F = are that (1) we have an inner product defined on all (p+ l)-forms 
F, and (2) that this inner product is positive definite. In Minkowski spacetime, the Hodge 
inner product [FA -kG is neither defined for all [p + l)-forms, nor positive definite, so even 
compactly supported p-form electromagnetic fields in the vacuum need not vanish. 

There is one way we might try to weasel out of having such a trivial classical theory in cases 
where conditions (1) and (2) are met. Namely, we can take a pre- gauge theory perspective on 
electromagnetism, and take F as the fundamental field of classical electromagnetism, rather 
than A. If the electromagnetic field F is not necessarily the differential of any gauge potential, 
then dF = is not automatically satisfied, so we must include it explicitly in the p-form 
Maxwell equations: 

dF = {} and d'^F = 0. 

Wc will see that these equations do have solutions other than F = 0, but only on spacctimcs 
with a certain topological properties — roughly, spacetimes with 'holes' of the type which 
can be enclosed by a surface of dimension p + \. What we will find is that classical p-form 
electromagnetism in Riemannian signature, in cases where we have an inner product defined 
on all {p + 1) forms, is what we might call a topological classical field theory: its space 
of solutions is determined up to canonical isomorphism by the topology of spacetime! This, in 
fact, is true whether we take ^ or F as the basic field, though the two perspectives disagree 
on what aspect of the topology determines the space of solutions. 

In the next section, we quantize M p-form electromagnetism and, in parallel, develop enough 
topology to state a criterion for the convergence of the path integral. As this topology — 
cohomology — is precisely what we need to fully understand the comments above about the 
classical theory being topological, we return to the p-form Maxwell equations in Section 6.2. 

6 Cohomology and Path Integrals 

6.1 A Criterion for Path Integral Convergence 

Consider some small portion of a complex of R-valued cochains 

dn—l da 

ci-^ CI — ^ c«+i 
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Cohomology classifies topological spaces by comparing two subspaces of C^: 

Z'i := kerdy — the space of g-cocycles 

B'l := ran dq^i — the space of g-coboundaries 

For the cochain complex of any spacetime lattice we have, by the defining equation dd = 0, 

C Zi. 

That is, every g-coboundary is a g-cocycle. Whether the converse of this statement is true de- 
pends on the particular topology of the spacetime lattice. If every g-cocycle is a g-coboundary, 
so that B*^ and Z*^ are equal, we say that the cochain complex is exact at C. In topologically 
interesting spacetimes (or regions of spacetime) , exactness may fail, and we measure the failure 
of exactness by taking the quotient space 

HI := ZyBi, 

called the gth cohomology with real coefRcients. 

To sec the ramifications of cohomology in lattice electromagnetism, let us reexamine what 
caused the path integral 

Z = [ e-^^^^VA 
Jcp 

to diverge in Section 2. In gauge theory, the action S{A) is invariant under gauge trans- 
formations. In our lattice version of p-form electromagnetism with gauge group M, gauge 
transformations take the form 

At-^ A + d(p, 

where ip is an arbitrary (p — l)-cochain. In other words, two p-connections A, A' G are 
gauge equivalent if they differ by a p-coboundary d(p G B^. In Section 2 we saw that the 
existence of gauge freedom implied path integrals such as the one above must diverge. Indeed, 
if ^ = BP is nontrivial, then we can write 

CP = BP(B (BP)^ 

and express any p-connection A uniquely as a sum 

A A , A± -.u ^0 e BP 
^ = ^o + ^^with ^!l^^^p^± 

This lets us break the path integral up as 

Z= [ [ e-^(^o+^")dA^dAo. 

J BP 7(BP)-L 

But since the action is gauge invariant, -S'(^o + A^) = S{A^), and so this becomes 



Z= / e-^'^^^UA^dAo 

J BP J{BP)-^ 

= [ dAo [ e-^^^^UA^ = oo, 

J BP J{BP)-^ 
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since for each fixed value of A-^, the integral over diverges. 

To eliminate these divergences it is therefore necessary that we eliminate any gauge freedom 
from our space of p-connections. The standard procedure is to pass to the quotient space 

A CP 

— = — = p-connections modulo gauge transformations, 

effectively declaring gauge equivalent p-connections to be equal. 

In many cases, integrating over C^/B^ instead of just does indeed make the path integral 
finite. But this does not always work! Whether it works or not depends on a topological 
condition — the cohomology of the lattice. While factoring out gauge freedom is necessary, it 
is only sufficient if the pth cohomology, HP is trivial. Indeed, if 

ZP 

BP ^ 

then since 

A _CP ^CP ZP 
'g ~ ^ ^ 'zp ® BP 

we have 



/ e-'^^UA= j [ e-'^^UA,dA2 

JA/g Jzp/BP Jc'P/ZP 

= I dA2 I e-^^^^UAi = GO, 



IZP/BP JCp/Zp 

since S{A) really only depends on dA = d{Ai + A2) = dAi. 

Pondering more carefully what caused the divergence, we realize that to eliminate all 
divergences from M electromagnetism we must use the quotient space 

A CP 

—r- = -=- = p-connections modulo flat p-connections 
Ao ZP ^ ^ 

For any lattice with only finitely many p-cells, integrating over this space always gives us a 
convergent path integral. The point is that, when HP is nontrivial, the quadratic form S is 
still degenerate, even after factoring out all of the gauge freedom. However, when we factor 
out not just gauge transformations, or in other words 'pure gauge' j3-connections, but all flat 
p-connections, S becomes a positive definite quadratic form, so that J e~^ is finite. Visually, 
on the space of p-connections modulo flat p-connections, every cross section of the graph of 
e~^ looks just how we expect a well-behaved Gaussian to look: 
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in contrast with the analogous picture from Section 2 



Let us consider two examples: one in which modding out by gauge transformations is 
sufficient to render the path integral finite, and one in which it is not. First, let us finally fix 
the problem we found in Section 2, namely the divergence of the path integral for 




over the space of connections. We write the connection as 

A = A-ye" + A^f? + A'if? = Aie^ 
where {e*} is the basis of 1-cochains dual to the basis {cj} of edges. A gauge transformation 

acts on the connection hy A A + del) where 

# = ((^2-</'i)(e^ + e2 + e3). 

Since the number (/>2 — (/>i € M may be chosen arbitrarily, the space of gauge transformations 
is just the 1-dimensional subspace 

:= ran do = span{e^ + + e^} C 

Similarly, the curvature F is given by 

F:=dA= {Ai - A2)p^ - {As - 

so the space of flat connections is 

^0 = Z'^ := kerdi = {A G Ci\Ai = A2 = A3} 
= span{e^ + + e^} = B^. 

Thus, due to the simple topology of the model, the flat connections in this case are precisely 
the pure gauge connections. In other words, the first cohomology is trivial: 

= Z^/B^ = {0}. 

Now let us do the path integral over the space C^/B^ = C'^ jZ^. To do this we make use 
of the canonical isomorphism 

^1 V ; 
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given by associating with each gauge equivalence equivalence class [A] G C^/B^ its unique 
member which is also an element of {B^)^. We first define a new orthonormal basis 



„3' 



1 

1 



y + e2 - 2e^) 



of chosen so that {e^ } is a basis for B^ and {e^', e^'} is a basis for (B^)-^. 
We can write the connection A in terms of the new dual basis: 



A = Ajie^' 



Then the action becomes 



^1+^2 + ^3 1' , A 



V3 



1—^2 2' 



V2 



A1 + A2- 2A 



3„3' 



^/6 



1 f{Ai-A2f , {A2-Asf 



+ 



Vo 



— [(A^-A2f + {A2-Asf] 



2e'^ViV2 
1 

2e^ViV2 
1 

2e'^ViV2 



V2{^A2<f + ^l(-^^2' + y|^3' 
[^2' ^3'] 



2V2 + ^Vx ^Fi 



^2' 
Ay 



Now integrating over AjQ is a snap! Since the 2x2 matrix in the above expression is non- 
singular, the formula we tried to use for this same problem in Section 2 now gives a finite 
result: 

jgf y det(Q) V 3 

As a second example, consider again the case of ordinary 1-form electromagnetism in two 
dimensions, but this time in a spacetime with the topology of a torus: 
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Note in the diagram that objects with the same label are identified. Let us first calculate the 
cohomology of this lattice. The cochain complex 



c° — > — ^ 



has differentials represented by the matrices 





-1 1 



1 -1 



10-10 
-10 10 



relative to the bases {v^,v^,v^} of C°, {e^, e^, e^, e^} of C^, and of C^. We are 

particularly interested in the first cohomology, since it relates fiat connections to gauge trans- 
formations. Some basic linear algebra shows that 

:= ker di = span{e^ + e^, e^ e^} ^ M 
:= ran do = span{e^ - e^} ^ M 



and therefore 



B^ M 



Physically, this says there are flat connections A on the torus which arc not gauge equivalent 
to the trivial connection, and that there are two degrees of freedom for such connections. In 
electromagnetism we also interpret = by saying that there are two Bohm-Aharonov 
modes [3, 29]. 

For the sake of characterizing the topology more fully, we should calculate the other coho- 
mology groups as well. We find 



and 



span{u^ — v^} 
span{p^ — p^} 



{0} 



For doing path integrals on this lattice, it is convenient to take the following orthonormal basis 
for C^: 



1' 



e = 



1 

7^ 



V2 



„3' 



V2 



V2 
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since then {e^ } is a basis for (Z^)-^ and {e^', e^', e^'} is a basis for (B^)-^. In terms of this 
new basis we can write the connection as 

A - A J - - ^3) y . {Al + ^3) ^2' , {A2±A^.' , {A2 - ^4) „4' 

\/2 \/2 \/2 

For simplicity, let us assume that the area of each plaquette is 1, so the action becomes 

SiA) = ^{{A, - Asf + (As - A^f) = ^Af,. 

Now just as in the previous example, integrating over connections mod gauge transformations, 
i.e. C^/B^ is the same as integrating over {B^)^ = spanje-*^ }, relative to which basis 

A has components {Ai' , A2' , Ay ). But S{A) is obviously a degenerate quadratic form in the 
variables Ai/, A2', A3/, so the path integral diverges: 

f e-^(^)PA = f e-^^^''/^'dArdA2'dAy = 00. 

Similarly, integrating over connections mod flat connections, i.e. C^/Z^ is the same as 
integrating over [Z^)-^. But S{A) is obviously nondcgcncratc in the variable Aii, which is the 
coordinate of A relative to the basis {e^ } of (Z^)^. In fact, in this case we get: 



Jc^/z^ Jm V 2 



So we see that to regularize the path integral for the torus, in M electromagnetism we 
must kill off both of the Bohm-Aharonov modes, in addition to factoring out gauge freedom. 
We have seen that the criterion in M electromagnetism for path integrals over the physical 
configuration space A/G to converge is that the first cohomology be trivial: 

\ converges for M electromagnetism / 

More generally, p-form electromagnetism has what we might call the p-form Bohm- 
Aharonov effect, but whereas the ordinary Bohm-Ahaxonov effect applies to regions of space- 
time with nontrivial first cohomology, the p-form version depends on the pth cohomology. We 
have also seen that getting path integrals in M p-form electromagnetism to converge requires 
killing off all 'p-form Bohm-Aharonov modes': 

\ converges for M p-form electromagnetism I 

In fact, when is trivial, we can sometimes simplify the process of taking path integrals 
further. To see this, note that the homomorphism dp : C^+^ induces an isomorphism 

ZP 
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In the case where HP is trivial, C^/ZP = C^/BP, which is the physical configuration space. If 
in addition dp is onto, so that B^'^^ = C^'^^, we then have an isomorphism 

g ZP ~ 

In other words, when there are no Bohm-Aharonov modes, and dp is surjective, we can just 

as well do path integrals over the space of curvatures F, rather than over the space of p- 
connections. In our first example above, is trivial and every possible curvature is d of some 
connection, so we could calculate the partition function more easily as: 

Z= f e-^(^''+^^'hFidF2 = 27re2 

The value of the partition function is different in this case, but this does not matter since Z 
is just a normalization factor — what matters is that expected values of observables are the 
same. 

For p-form electromagnetism with gauge group M, we thus have several choices for what 
space to integrate over when we do path integrals: 

1. A, the space of p-connections — This is the naive approach of Section 2, and works 
almost never: only when there are no nontrivial flat connections. 

2. A/ G, the space of p-connections mod gauge transformations — This works whenever there 
are no p-form Bohm-Aharonov modes, i.e. when spacetime has trivial pth cohomology. 

3. A/Aq, the space of p-connections mod flat p-connections. This always gives convergent 
path integrals, and is the same as A/Q when the pth cohomology is trivial. However, in 
the case of nontrivial cohomology, it means ignoring the p-form Bohm-Aharonov modes. 

4. := C^+^, the space of curvatures of p-connections. In the case where the pth cohomol- 
ogy is trivial and d: CP ^ (7^+^ is onto, this is equivalent to the previous two options, 
and generally far easier to calculate with. 

5. J-Q := BP^^, the space of curvatures F, subject to the constraint dF = 0. This is like the 
previous option, but works even when d: CP ^ (jp+i jg j^q^ onto. The practical tradeoff 
is that the constraint makes calculating integrals less straightforward. 

The major shortcoming in the case of gauge group M is that none of these options allows us 
to take any p-form Bohm-Aharonov modes into account and get convergent path integrals. 
We see in Section 7 that this problem is resolved when we switch the gauge group to U{1). 
We remark that options 4 and 5, while sometimes convenient in the free held setting, are less 
viable when we add matter, since matter fields typically couple to the A field, not F. 

Let us look at some examples. We have seen that the torus has nontrivial cohomology. In 
fact, for purposes of computing cohomology, we could have trimmed down our example to a 
torus with only one vertex, two edges, and one plaquette: 

Vl ei Vl 



'62 



Vl ei Vl 
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Since cohomology is a topological invariant — topologically equivalent spaces have isomorphic 
cohomology — we find again the cohomology groups 



2 



In particular, H tells us there are two Bohm-Aharonov modes, just as there were for our two- 
plaquette torus above, and just as there would be for any other spacetime with the topology 
of a 2-torus. 

As another example, consider a discrete spacetime with the topology of an n-sphere, S'^. 
When the gauge group is M this has cohomology groups 



Thus, for example, there are no Bohm-Aharonov modes for ordinary electromagnetism on the 
2-sphere, since H^{S'^) = 0. Yet for the sphere the differential di is not surjective. Thus in 
this case we may do path integrals over the physical configuration space A/G, but not over 
the entire space of curvatures, C^. By contrast, for a Kalb-Ramond field on the 2-sphere, 
there are 2-form Bohm-Aharonov modes , since H'^{S^) = M. In this case, R-electromagnetism 
gives divergent path integrals unless we kill these modes. But what about the other nontrivial 
cohomology, H^(S^)7 This evidently corresponds to a Bohm-Aharonov effect for '0-form 
electromagnetism'. We consider this in Section 8. 

6.2 Hodge's Theorem, Topological Classical Field Theory 

Having quantized lattice p-form electromagnetism, let us step back and take a closer look at 
the classical theory. Recall from Section 5 that the vacuum discrete p-form Maxwell equations: 



have two possible interpretations: 

1. Considering the gauge potential A as fundamental, we have seen the positive definiteness 
of the inner product on (p + l)-forms reduces the equations of motion to F = 0. From 
this perspective, then, solutions of the vacuum discrete p-form Maxwell equations are 
fundamentally p-forms A satisfying dA = 0, modulo exact p-forms. In the global lan- 
guage, solutions are precisely elements of the space of fiat p-connections modulo gauge 
transformations. But this is just 



the pth cohomology. 

2. Considering the electromagnetic field strength F as fundamental, solutions of the vacuum 
discrete p-form Maxwell equations are (p -|- l)-forms F satisfying 




d^F = and dF = 



g BP 



dF = and d^F = 0. 
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The second of these equations says precisely that F is orthogonal to every exact (p + l)- 
form, since the nondegeneracy of the inner product says d^F = if and only if {SF, X) = 
for every p-form X . or in other words {F, dX) = for every exact {p+ l)-form dX. We 
show below that positive definiteness of the inner product on p + 1 forms implies that the 
subspace of closed {p + l)-forms F which are also orthogonal to every exact (p+ l)-form 
dX, is isomorphic to 

the + l)st cohomology. 

So, from either viewpoint, the classical theory is 'topological'; but the two viewpoints disagree 
on which aspect of the topology determines the solutions! 

To prove the claim in item 2 above, it is helpful to develop an alternate viewpoint on 
cohomology — a viewpoint which is useful whenever we have a cochain complex of (possibly 
indefinite) inner product spaces. When we have such a complex C, each of the differentials 
dq-.Ci^ has an adjoint 4 : Ci+^ C?: 

1 d d 

ci-^ — ^ c« - — ^ C«+i 

dt (it 

The first important observation is that dd = together with nondegeneracy of the inner 
product implies d^d^ = 0, so the maps d^ also form a complex. We say w G is coclosed 
if d^LO = or coexact if u = d^a for some a G C^"^^. Just as we noted in the special case 
of F above, a coclosed g-form is precisely a g-form which is orthogonal to all exact g-forms. 
Similarly, a closed g-form is precisely a ^-form which is orthogonal to all coexact ^-forms, by 
the same argument. That is: 

kerdq = (ran o?^)"*" kerdg_^ = (ran dq-i)'^. 

We shall combine these to prove the Hodge decomposition^ which says we can write as 
the orthogonal direct sum 

= ran dq-i © ker Aq © ran d^. 

Here 

Aq := d\dq + dq^id\_^ 
is the Laplacian on C"^, and a cochain in its kernel is said to be harmonic. 

Proof of the Hodge Decomposition: The remarks immediately preceding the statement 
of the Hodge Decomposition give us two separate direct sum decompositions of C^: 

= ker dq © ran d| = ker ^ © ran dq-i. 

^There arc more sophisticated versions of the Hodge decomposition, such as the usual version on the dcRham 
complex, and the 'Kodaira decomposition' for complexes of operators that are only densely defined [7]. We supply 
the proof of this simple version here because it is illuminating. Proofs of fancier versions are the same in essential 
structure, though some nontrivial analysis is involoved. 
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The first of these says we can write a given u £ uniquely as a; = u;o + with loq G ker dq. 
We then use the second direct sum to write ujo = ujh + da with u>h € ker But in fact, this 
Ufi = loq — da is also in the kernel of dq since loq and da both are. We thus have 

(jj = da + (jOh + d^ P with u>h G (ker d]^_^) fl (ker dq) 

Moreover, this decomposition is unique, since iOh+d"^ (3 and da are uniquely determined by C"^ = 
kerd^_j^ © ran dq-i, while ujh + da and d^ P are uniquely determined by C = ker © ran dq. 
The sum is thus direct: 

C = ran d © (ker d^ D ker d) © ran d^ , 

and orthogonal since uj^, da, and are mutually orthogonal. 
Now suppose Lo is harmonic. Then 

= {u;,{dU + dS)u) 
= {d^u,d^uj) + {du,duj). 

Up to this point, we have not used that the inner product is positive definite. Using this fact 
now, we conclude that 

du = and d^uj = 0. 

Conversely, if a; solves these equations, then obviously {d'^d + dd^)^ = 0, so u is harmonic. We 
thus see that the space of harmonic g-forms is the space of closed and coclosed g-forms, and 
the decomposition is established. ■ 

Recall that the free p-form Maxwell's equations for F simply say that F is closed and 
coclosed. The above proof thus shows that when we have a positive definite inner product 
defined on all (p+l)-forms, solutions of Maxwell's equations are precisely the harmonic (p+l)- 
forms. In fact the Hodge decomposition lets us prove a version of Hodge's Theorem, which 
says the space of harmonic g- forms is canonically isomorphic to the cohomology H'^. To see 
this, write a given g-form uj as 

oj = da + coh + d) P 

with u)h harmonic. To understand cohomology in this context, we need to understand what 
it means for u to be closed. But we have observed that a closed form is precisely a form 
orthogonal to every coexact form d^- Since the Hodge decomposition of u is orthogonal, this 
means uj is closed precisely when its coexact part d^ P vanishes. Thus, noting that ran dq-i is 
just B'i, the Hodge decomposition 

= Bi® ker © ran d^ 

implies 

Zi = B'i®kerAq 

and hence 

HI := — ^kerAo. 
Bi * 

That is, the qth cohomology H'^ is canonically isomorphic to the subspace of harmonic q'-forms. 
Identifying these two spaces, we can therefore write: 

CI = dC^-^ © F« © d^Ci+^ 
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as an alternate version of the Hodge decomposition of C"^. 

In the case we are most interested in, the case q = p + 1, we thus see that the solutions of 
the p-form Maxwell equations 

dF = and Sf = 

are 

ker d n ker d^ = ker Ap+i ^ HP^^. 

That is, when we consider F as the fundamental classical field, solutions of the equations of 
motion are specified by the {p + l)st cohomology H'p^^: 

dF = and d^F = ^ F e HP+^. 

7 U{1) p-Form Electromagnet ism and the p-form 
Bohm-Aharonov Effect 

7.1 Gauge Groups for p-form Electromagnetism: M vs. U{1) 

We have seen that eliminating all of the divergences from path integrals in R p-form elec- 
tromagnetism requires that we factor out not only all of the gauge freedom, but also the 
Bohm-Aharonov modes. This is undesirable: the Bohm-Aharonov effect is an empirical fact 
in ordinary 1-form electromagnetism. Neither should we rule out the p-form Bohm-Aharonov 
effect in our calculations. The ultimate solution to the problem of divergent path integrals is 
to switch our gauge group from M to U{1). We shall see that U{1) electromagnetism cures 
all divergences, even in cases where the pth cohomology is nontrivial! Better yet, C/(l) path 
integrals converge even before we factor out the gauge freedom. The reason for this is simple. 
The space of U{1) p-connections is 

A{U{1)) = U{1)^^ 

where Xp is the set of p-cells in the lattice. But this is just a product of circles — a torus! 
Since a torus is compact, 

converges for any continuous function / of A. Of course, for / to be a physical observable, we 
still want / to be gauge invariant: f{A) = f{A + d(j)). 

Indeed, the switch to U{1) promises to be such an improvement that the reader may well 
wonder why we have bothered with the M case. As it turns out, the U{\) theory can be 
obtained most easily from the M theory. In particular, without having developed the M theory, 
it would be difficult to guess the best action to use in evaluating U{1) path integrals. The 
reason is that since the curvature lies not in the vector space M^^'+i but in the mere group 
we do not have the same analogy between cochains and and differential forms as we 
had in the M case. In particular, we cannot rely directly on an inner product of cochains for 
the action as we did in the M case in Section 4. What we will show is that we can turn p-form 
electromagnetism with gauge group R into p-form electromagnetism with gauge group U{1), 
essentially by "wrapping the real line around the circle." 
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The key to understanding the relationship between the abeUan gauge groups M and U (1) 
is the following 'short exact sequence' of homomorphisms: 



^ Z > R ^ U{1) ^ 

where the second map is the usual inclusion and the third sends x G R to e^'^*^ G U{1). Here 
'exactness' means the same as for the chain and cochain complexes: that the range of one map 
equals the kernel of the next. In particular, this says we map R to U{1) by winding the line 
around the circle once every 2tt, and that the kernel of this map is precisely Z. 

In comparing different gauge groups, it is convenient to write the group of /c-form fields 
not simply as C'^ but more explicitly as 

C\G) = hom{Ck,G), 

i.e. homomorphisms from the group of fc-chains to the gauge group G. 
Now given any R p-connection 

we get aU{l) p-connection 

i:= e^^^^: Cp^U{l) 
by composition with the homomorphism R — U{1). We thus get a map 

hom(Cp,R) ^ hom(Cp, U{1)) 

and the fact that the above sequence of abelian groups is exact implies that the kernel of this 
map is hom(Cp,Z). Moreover, since Cp is free, every U{1) p-connection comes from an R 
p-connection in precisely this way. Given a p-connection A: Cp U{1), map each generator 
c in the basis Xp of Cp to some A{c) G R with eicp{2mA{c)) = A{c). This defines an R 
p-connection A which 'lifts' A: 

Cp 

y" \a 

R^^f/(l) -0 

That is, we have the identity: 

A = A. 

In fact, each of the sequences 

^ hom(Cjk, Z) ^ hom(Cjk, R) ^ hom(Cfe, U{1)) ^ 
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is exact, thus giving a short exact sequence of chain maps: 

^ hom(C„, Z) ^ hom(C„, M) ^ hom(C„, f/(l)) 



^ hom(Ci, ] 
d 

hom(Co,] 

which we may write more succinctly as 

> hom(C,, Z) ^ hom(C., 1 



0- 



■hom(Ci,Z) 
d 

■hom(Co,Z) 



■hom(Ci, [/(!)) 

d 

hom(Co,f/(l)) 



hom(C.,i7(l)) 



-^0 . 



Lattice p-form electromagnetism with abehan gauge group G is all about the cochain complex 
hom(C,,G). We may thus express the physical content of the above exact sequence of chain 

maps by saying there is a projection from the theory with gauge group M to the theory with 
gauge group U (1), as we wanted, and that the kernel of this projection is a theory with gauge 
group Z. Metaphorically: 



/ lattice p-form \ 
electromagnetism 
with gauge group 
Z 



/ lattice p-form \ 
electromagnetism 
with gauge group 



V 



/ lattice p-form \ 
electromagnetism 
with gauge group 

m 



-^0 



is exact. 

One critical implication of the exact sequence of cochain complexes is that two M p- 
connections A, A' G hom(Cp,M) which are gauge equivalent, say 

A-A' = d(t> ^ € hom(Cp_i,M), 

project down to ^7(1) p-connections A, A' G hom(Cp, ^7(1)) which are gauge equivalent: 

A-A' = d4> 4>ehom{Cp-i,U{l)). 

Briefly, M gauge equivalence implies U{1) gauge equivalence. 

But in some cases, there is a genuine difference between the two choices of gauge group. A 
simple example is a space which is sometimes called an 'm-fold dunce cap', DC„i. This can 
be constructed from one vertex v, one edge e, and one plaquette p which is sewn along its 
boundary around the looped edge in such a way that it wraps around m times: 
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The dunce cap has cohomology groups 

H°{DCm,'R)='!^ H^{DCm,'^) = ii"2(DC^,R) ^0 

for G = M, but 

H''{DC,n,U{l)) = U{l) H\DC,n,U{l))=Z/m H\DC^,U{1)) = 

for G = U{1). The significant difference is in the first cohomology. The fact that is trivial 
in the M case says every R-valued 1-cocycle is a 1-coboundary. Said another way, all flat M 
connections on DCm are gauge equivalent: 

dA = dA' = ^ d{A-A') = =^ A-A' = dcj). 

In fact, this result is not surprising given our particular construction of DCm — there is only 
one flat M connection, the zero connection. But in the U{1) case, is a cyclic group of order 
m. It is not hard to see why. There are m flat U{1) connections on DCm, given by assigning 
to the edge e any one of the m distinct mth roots of unity in ^7(1). 



7.2 Gaussian Integrals on a Torus 

In the U{1) version of lattice electromagnetism, the holonomy of the connection takes values 
on the circle U{1). To do path integrals, we need the analog of a Gaussian on a circle, which we 
obtain from the ordinary Gaussian by wrapping the real line around the circle. In particular, 
since 



is a probability measure on 



27r(7 



(1) 



is a probability measure on = U{1) which we shall refer to as a circular Gaussian. We 
can use this measure to calculate the expected value of a function / : U{1) — >Mof the random 
variable 6 as follows: 

/•27r -(e+2n7r)V2a2 >/, 

/we ; 

f2(n+l)7r 



1 pzin+ijir 



2na 



f{9)e-'"/'"'"de 



What this calculation shows is that integrating a function /: U{1) M against the measure 
in (1) is the same as extending / to a periodic function / : M — > R and then integrating this 
against the usual Gaussian measure. 



38 



Since a periodic function on M (or equivalently, a function on the circle) can be expanded 
in a Fourier series, we should work out the expected values for each element of the Fourier 
basis {e'"*^}. For any m G Z we have 



1 

y2Tra J-oo 
V^ira J- 



e 



2'Ka J-oo 

t / g-(6'2-2im(T20-m2(T4)/2CT2g-m2 0-2/2^^ 

2'Ka J-oo 

—m^a^ 12 rco 
J-oo 



'l-KG J -oo 
g-m2aV2 



,-m2aV2 



27rC7 



Notice that this expected value is real, even though e*"^^ is complex. The imaginary part is 
zero on average, since the imaginary parts sin(m^) are odd about the circle's origin. 

In fact what wc arc calling circular Gaussians are actually famous special functions in 
their own right, though they are usually presented in a slightly different form. Using the 
Fourier coefficient ^e*™^) above, one can expand the circular Gaussian itself in the basis {e*™^}, 
resulting in the identity: 

g-(e+2n7r)2/2a2 = J_ ^ ^-n^ay2^ine_ 



The form on the right side of this equality makes it easier to recognize that our circular 
Gaussian with 'standard deviation' is really: 

2iT \2tt' 27r J 

where 

neZ 

is the theta function [24]. 

When it is convenient, we continue to refer to this theta function loosely as a 'circular 
Gaussian', and write it though this is not technically the exponential of 

anything. We engage in this gross abuse of notation because it is suggestive of how one actually 
calculates integrals using the corresponding 'theta measure' — by expanding the circle to a 
line and integrating an ordinary Gaussian. 

The higher dimensional version of the above calculations follow the same form. Suppose 
we have a Gaussian 
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on R^. Following the prescription we used for the circular Gaussian above, we wrap each of 
the N coordinates around a circle: 



e 2 



to obtain an analog of a Gaussian on a torus, what we might call a toroidal Gaussian. Then: 



/det(Q^ r imHu^-\Q'ie,e,.n 
(27r)^ 7m- 

/det(Q*J)^_iQ-i^fc^i /■ ^-hQ'H0i0j-2iQJ^^m'^ei-Q-^'QT^^m''m')^Q 



(27r) 



fdet(Q^i,^__ 



Using this and proceeding as in the 1-dimensional case, we obtain the identity 



,2 



which reduces to (2) in the case N = 1, Q = 1/a . The right side is again a theta function, 
this time its multivariable generalization: 



(27r)^ V^tt' 27r 
Here 

'd{z,n) := ^7rin-Qn+27vin-z 

where z G C^, and is a symmetric N x N complex matrix with positive definite complex 
part. [24] Again, wc continue to write this theta function as "e~2^''3^" and think of it as a 
'toroidal Gaussian' when convenient, since this is evocative of the analogy to Gaussians on 



7.3 U{1) Path Integrals 

Recall that for gauge group M, our action for discrete p-form electromagnetism is given by 

2>'^) = 2i-- 



SiA) = ^{F,F) = ^h^^F,Fj 
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where h''^ is the matrix of the inner product relative to the basis of (p + l)-cochains consisting 
of dual {p+ l)-cells. As discussed in the previous subsection, making the transition from R to 
U{1) involves replacing the ordinary Gaussian exp(— 5") in the real variables Fi by a 'toroidal 
Gaussian' in the C/(l)-valued variables Ff. 



-S{A) 



e 2e 



^ h"^ [Fi -2niTv){Fj -2nj'iv) 



e 2e 



^h*^A{dxi)-2niTr){A{dxj)-2njTr) 



(4) 



where A'' = |Ap+i| is the number of {p + l)-cells in the spacetime, which we now assume to 
be finite. Note that we do not define the action S but only the analog of its exponential e~^, 
which is all we need for doing path integrals. 

For the sake of comparison, it is perhaps best at this point to return to our old two-plaquette 
spacetime from Section 2: 




The corresponding chain complex is still as constructed in Section 3 (and later used in Section 
6.1), but our gauge field A is now a C/(l)-valued 1-cochain: 

Aehom{Ci,U{l)). 

Thinking of C/(l) as the real numbers mod 27r, if 

[Ai] e M/27rZ ^ C7(l) 

is the value of A on the edge Cj, then using the action from Section 4.2 in equation (4) we 
obtain 

g-5(A) ^ ^_(A2-Ai-2ni7r)V2e2Vi-(A3-A2-2n27r)2/2e2V2 
neZ2 

where A^ is any representative of the class [Ai]. Conveniently, the fact that the matrix of the 
inner product is diagonal in this case lets us split up into one term for each plaquette: 

g-5(A) ^ g-S(A|pJg-S(A|p2) 



where 



g-5(A|pJ ^ ^ ^_(A2-Ai-2n7r)V2e2Vi 

g-S(A|p2) ^ ^ g-(A3-A2-2n7r)2/2e2y2 
neZ 
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Note the correspondence: since the path integral in the gauge group M case involves a product 
of Gaussians, one for the holonomy around the boundary of each plaquette, the path integral 
for the U{1) case involves a product of theta functions — the 'circular Gaussians' we defined 
above. 

We are ready to calculate the partition function for our example: 

f2n f2TT f2TT 



Jo Jo Jo mez 
Jo Jo ■^-^ 

Jo J -OO 

= (27re)2yTY^ 

The partition function is finite, in contrast to the M case. 

In the U{1) case there is actually a nice geometric interpretation of the coupling constant 
1/e^. Looking back carefully at the relationship between Gaussians and theta functions, wc sec 
that the square of the electric charge may be reinterpreted as the radius of the circle U{1). If 
we take this radius to be effectively infinite, so that the circle becomes a line, then we are back 
in the realm of M electromagnetism with unit electric charge. This helps explain the factor 
(27re)^ in our result for the partition function above: as the electric charge becomes infinite, 
the partition function effectively diverges, as we verified is Section 2 that it does in the M case. 

It is worth noting that the split of e""^ into the product ]^e~'^*^^») of plaquette actions 
happens more generally whenever the inner product is diagonal, and in particular in the case 
of p-form electromagnetism in p + 1 dimensions. The partition function for U{1) p-form lattice 
electromagnetism in the special case of p + 1 dimensions is thus 



= /" TT g-(A(ac)-2n7r)V2e2vol(c)p^ 



8 The Free Scalar Field and 0-Form Electromagnetism 

The alert reader may have noticed that the formalism for lattice p-form electromagnetism 
developed in the previous sections does not seem to apply equally to the case p = 0. In 
particular, there seems to be no notion of gauge invariance in the 0-form case. By extending 
our notion of n-complex we now correct this discrepancy. Lest the reader think this an esoteric 
digression of purely mathematical interest, wc note that 0-form electromagnetism is really just 
scalar field theory! Indeed, the 0-form Maxwell equation 

Sdcf) = 
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really just says, in the Riemannian signature we have been using, 

= 0, 

or in the more familiar Lorentzian case: 

□^ = 0. 

The field (f) in free scalar field theory — the 'gauge field', though the theory is too simple to 
be considered a 'gauge theory' by most — is just a G-valued 0-cochain in the lattice context: 

(f) G hom(Co,G), 

i.e. it assigns to each vertex an element of the gauge group G = M or But this is just what 

we called a 'gauge transformation' in electromagnetism. To figure out 0-form electromagnetism, 
let us therefore think again of our motivating example — ordinary 1-form electromagnetism, 
taking the gauge group to be M again for simplicity. 

One key point about gauge transformations in electromagnetism is that different connec- 
tions can induce the same curvature. But further thought along these lines reveals that this 
phenomenon is not at all unique to connections. Indeed, if we move one notch down the chain, 
different scalar potentials may give rise to the same connection. Texts on electromagnetic 
theory usually summarize this observation by saying the potential is only defined up to an 
additive constant. In fact, adding a constant to the potential is like "gauge transforming the 
potential." Better yet, it is like making a "mcta-gaugc transformation between gauge trans- 
formations." People don't usually talk about a choice of origin of the potential in these terms, 
perhaps because it is too simple, but this really is a primitive example of the 'ghosts of ghosts' 
one might read about in string theory papers. 

To make this idea precise, we use the relationship between the connection A and potential 
(or scalar field) ^ as a model for that between potential and what we might call prepotential. 
This seems to demand we add a new term, some C~^, to our cochain complex, along with a 
new differential d_i : C°. Then a meta-gauge transformation of a scalar potential (f) by 

a meta-potential u G should look like: 

(p 1-^ (p + d-iu. 

What should and d-i look like? Usually we only think of adding a constant to the 
potential, but that is really only because of our peculiar fondness of connected spacetimes. 
In general, we can add a local constant to (f) without affecting dcj), so can have as many 
degrees of freedom as there are connected components in the lattice. In fact, we might as well 
take C-i to be the free abelian group generated by the set of connected components of the 
lattice, and let C"^ = hom(C_i,M) be its dual. 

Now to make i— (j) + d-iu act as an adjustment of the origin of ^ in each connected 
component, we simply define d^iu: Cq ^ K to be the map which assigns to each vertex (i.e. 
each basis element) in Co the value of u: C_i — > M on the connected component in which that 
vertex lies. That is, given a vertex v G Co, if we denote its connected component by dv, then 
d-\u{v) := u{dv). This, of course, defines a new boundary map d-i : Co — C_i — it is just 
the adjoint of cZ_i. 



43 



Putting all of this together, we arrive at a new 'augmented' chain complex and cochain 
complex 

„ 9-1 ^ <9o ^ (9l ^ 02 dn ^ 

U < U_i < C-o < L/l < • • • L/rt 

What we have just seen is that the natural meaning of the notion 'boundary of a vertex' 
is just the connected component in which the vertex lies. It is interesting to note that at this 
level the old motto "the boundary of a boundary is zero" just says that the two cndpoints 
of an edge lie in the same connected component! As a final remark, note that we cannot 
continue extending the complex further left. Since each connected component contains (i.e. 
is the boundary of) some vertex, Oq is onto, so dd = implies the range of 9_i must be 
trivial. Phrased from the dual perspective, there can be no nontrivial gauge transformations 
in '(— l)-form electromagnetism'. 



9 Discretization Independence in p + 1 Dimensions 

Recall from Section 4.2 that when the gauge group is M and spacetime is (p + l)-dimensional, 
the action in lattice p-form electromagnetism in p + 1 dimensions takes a particularly simple 
form. Namely, it just involves the sum of = {dA)^ over all (p + l)-cells, weighted by the 
volume of each cell. 

SiA) = ^ V ^ 
^ ' 2e2 ^ vol(x) 

In this case we have the beautiful result that the theory does not depend on the details of 
the discretization! In fact, there is a strong sense in which the quantum theory of p-form 
electromagnetism in p + 1 dimensions is almost topological: it's only nontopological degree 
of freedom is the total volume of spacetime. It is thus what we might call a volumetric 
quantum field theory. 

To begin seeing why the p-form theory in (p+l) dimensions is volumetric, consider splitting 
a (p + l)-cell of volume V in two by slicing through it with a new p-cell, leaving some fraction 
of the original volume on either side. Schematically: 




V 




Vi + V2 = V 

We wish to compare the curvature Fq on the original (p + l)-cell to the sum F = Fi + F2 in 
the finer discretization. According to the action above, the curvature on the original cell is a 
Gaussian random variable with probability measure 

HFo = — r^dFo 



44 



where we suppress the couphng constant l/2e^ in this section to avoid cluttering our calcula- 
tions unnecessarily. 

On the other hand, the joint probability measure for the random variables Fi and F2 is 

tiFi,F2 = 7== dFidF2. 

Making the change of variables 

F = Fi + F2 
G = Fi — F2 

we can rewrite this as a measure for F and G: 

^^■^ = WTO 

where the 2 in the denominator comes from the Jacobian 



d{Fi,F2) 



d{F,G) 



To get a measure for F, we now have only to integrate out the G-dependence from iJ:f,g- 
Collecting terms, using V = Vi + V2, and then completing the square in G, we get 

^iFG = ^=e ^^1^2^ ^ ' dFdG 

Integrating out the dependence on G then gives a normalized measure for F = Fi + F2: 

fiF = — dF, 



which is the same as ij,Fo- 

What this calculation shows is that we are free to rediscretize our spacetimc by knocking 
down or inserting new (j) + l)-cells, and that we will always get the same results, as long as we 
only ask questions that can be asked for either discretization. In the end, the only degree of 
freedom besides the topology should be the total volume of the spacetime. This result could 
be made more precise. To complete the argument, one would need to establish a complete set 
of 'moves' which allow passage from one discretization to any other, and then show that the 
theory is invariant under these moves, as we have done for the move described above. Since 
it is special to the case of p + 1 dimensions, however, we instead focus on a more general 
description which works in arbitrary dimensions. This is the subject of the next section. 



10 Chain Field Theory 

Electromagnetism is not a topological quantum field theory. In Section 9 we began to see that 
the special case of pure p-form electromagnetism in {p + 1) dimensions is almost topological 
— it requires only one nontopological datum: the volume of spacetime. But it is interesting 
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to see to what extent our more general theory of p-form electromagnetism in n-dimensional 
spacetime (including ordinary electromagnetism in 4 dimensions) can be given a TQFT-like 
description. If we pick an inner product on M-valucd {p + l)-cochains and define the action 
as in Section 4, then this goal can indeed be realized to a large extent. Before diving in, we 
warn the reader that this section requires a slight jump in mathematical sophistication. In 
particular, we need to assume an acquaintance with category theory, referring the reader to 
reference [20] for details not presented here. 

In topological quantum field theory, 'space' at any given time is a compact (n — 1)- 
dimensional manifold, whereas spacetime is an n-dimensional cobordism connecting two slices 
of space. This is easiest to visualize when n = 2, where space is necessarily just a union of 
circles, and a cobordism between one slice of space and another is a 2-manifold with boundary^ 
having these two slices of space as the two components of its boundary, like this: 




A TQFT assigns to each (n— l)-manifold S representing space, its Hilbert space of states Z{S), 
and to each cobordism M : S ^ S' between (n — l)-manifolds a linear operator Z(M) : Z(S) — ^ 
Z(S'), which we can think of as a 'time evolution' operator. 

More precisely, a TQFT is a symmetric monoidal functor [20] 

Z : nCob ^ Hilb 

from the category of smooth n-dimensional cobordisms to the category of Hilbert Spaces. 

Since we've described discrete spacetime in electromagnetism using chain complexes, we 
would like to define a chain field theory to be a symmetric monoidal functor 

Z : nChain Hilb 

where nChain is a category patterned after the category nCob, but with manifolds replaced 
by chain complexes. However, as we have already noted, electromagnetism is not merely 
topological. Consequently, the category nChain will have a richer structure than nCob. In 
particular, p-form electromagnetism as described in the previous sections relies on having an 
inner product on (p -|- l)-forms. 

Accordingly, an object S in the category nChain should be an (n — l)-complex (in the sense 
of Definition 1):^ 

■< Sq ■< Si ■< • • • ■< Sn-i 

^Morc properly, a cobordism is a diffcomorphism class of such 2-manifolds with boundary. 

''Here we begin a notational practice which we continue throughout this section. Since we will be dealing with 
multiple complexes, representing both space and spacetime, we denote, for example, Cfc(S', M) or Ck{S,U{l)) by the 
shorthand Sk rather than Ck- 
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equipped with an inner product on M-valued (p + l)-cochains: 

{ , )s- hom(S'p+i,M) X hom(S'p+i,M) M. 

To make our path integrals well defined, we demand that an object S have finitely many cells, 
so that, in particular, the dimension of Sp+i is finite. 

To define a morphism from S to S' Si nChain, let M be an n-complex, also with finitely 
many cells, whose chain complex M, is equipped with injective chain maps 



S, 



M, 



sending basis elements to basis elements. That is, for each k G {0, 1, . . . , n — 1} we have sub- 
group inclusions i: ^ and i' : S'^^ such that the following diagram is commutative. 

• • Sn-l 




■M 



n-l 



Mr, 



^ '-'n-l 

We also equip M with an inner product on M-valued {p + l)-cochains: 

( , )m- hom(Mp+i,R) X hom(Mp+i,M) R. 

and demand that M preserve the inner products from S and S' in the following sense. 
First, since Mp+i is free, we can think of SP~^^ = hom(S'p+i,M) as a subgroup of MP~^^ = 
hom(Mp_|_i, M) simply by extending any u G S^~^^ trivially to M^"*"^. That is, u becomes a 
member of M^'^^ by letting it be zero on MP'^^\S'P'^^ . We then demand that the inner product 
{—■>—)m is precisely {—,—)s when restricted to S^'^^ x S^'^^. 

A morphism in nChain will be defined as a certain equivalence class of the gadgets described 
in the previous paragraph. Suppose we have two such gadgets from S to S': 



•M, 



and 



s, m; ^ 



sL 



We say M and M' are equivalent, written M ~ M', if there exists an isomorphism of chain 
complexes 

(p: M, ^ 

preserving both the chosen bases of the and the inner product on M^, such that the 
diagram 

S, M, S', 



Ml 
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commutes. Finally we define a morphism 

[M]: S' 

in nChain, which we call a chain cobordism, to be the equivalence class under the relation 
~ of all gadgets S, — > M, < — S',. Following standard practice, we will generally drop the 
brackets from [M], writing simply M for our chain cobordism, blurring the distinction between 
the equivalence class and a representative. 

Now if M: 5 — ^ S" and N : S' ^ S" are chain cobordisms, we must describe a rule for 
composing them. Our mental picture of composing two cobordisms between manifolds as 
stacking them on top of each othersuggests composing M and N by stacking them to form the 
larger diagram 




■M 



n-1 



^ '-'n-1 
j 

f 

Of course, while this image is helpful conceptually, this diagram does not give us a new 
chain cobordism according to the definition given above. When we compose two cobordisms 
between manifolds, we don't just stack them on top of each other; we must glue their boundaries 
together in such a way that we get a new cobordism — in particular, a manifold with boundary! 
Similarly, when we compose our chain cobordisms M: S ^ S' and N: S' ^ S", we must get 
a new cobordism NM : S ^ S" 



S'J 



U (NM), 4 

by "gluing M and N together along S." Thinking of M and N as the physical lattices whence 
they derive, we want NM to have as n-cells all of the n-cells of both M and N. This leads us 
to define 

{NM)n := Mn © Nn. 

Likewise, NM should have all lower dimensional faces of both M and N, except that we 
identify cells in M and N if they come from the same cell in S under the gluing maps 



M, 



■N, . 



The natural way to accomplish this is to let 

Mk®Nk 



{NM)k :-- 



{i'{c)-j{c)\ceS',) 



k = 0,l,. 



,n 



1. 
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In words, each NM^ is the direct sum ® modulo relations that say c G C © 
and c. e Nk C © A?fe are equivalent if they are equal as elements of S'^. In particular, the 
diagram 

s:^— iv. (5) 

i' 

M, ^ {NM), 

commutes, where the unlabeled arrows are the obvious injections. Moreover this is the best 
way to make such a commutative square, in the sense that if K, is any other chain complex 
of abelian groups equipped with chain maps a and (3 such that 

i' a 

commutes, then there exists a unique chain map (NM), — making 



S', ^A^. 




commute. Category theorists express this 'universal property' of NM by saying the diagram 
(5) is a pushout square [20] . We mention this here because the pushout property streamlines 
the proofs of Theorems 1 and 2, below. Note that the {NM)k are still /ree, since the direct 
product of free abelian groups is free and identifying two basis elements of a free abelian group 
gives a free abelian group. 

We must define the inner product on (p + l)-forms on the composite NM. Since NM 
consists of all of the cells of N and M, there is a canonical way to do this. Given F,GE: 
hom((iVM)p+i,M), we define: 

{F, G)nm = {F\n, G\n)n + {F\m, G\m)m - {F\s' ,G\s')s' (6) 

Here the first two terms are just the inner products on M and iV, while the third term corrects 
for double counting on the common boundary S' of M and N . Specifically, if Xi and Xj are 
both {p + l)-cells in 5', then since S' is a subcomplex of both M and AT, each of the first 
terms in the inner product contains a term of the form (x*, x^). Subtracting {F\si, G\si)s' thus 
eliminates the extra {x^,x^) term. 

The reader already convinced, or else willing to accept on faith, that raChain as we have 
described it does indeed form a category may skip ahead to Theorem 2, where we prove the 
main result. For the more rigorously oriented, we collect some of the remaining technical 
details in the proof of the following theorem. 
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Theorem 1 nChain is a symmetric monoidal category. 



Proof: To show nChain is a category, we must show that composition of morphisms is well- 
defined (i.e. that it respects equivalence classes), that it is associative, and that identity mor- 
phisms exist. 

Suppose M ~ M' are two representatives of a chain cobordism from S to S', and N N' 
representatives of a cobordism from S' to S" . We must show that NM ~ N'M'. By the 
equivalence relation ~, there exist basis-preserving isomorphisms 



0: M, M', 

such that the following diagrams commute: 



S, 




■M, 

ML 




S', — ^ A^. S'^ 





Here is where pushouts come in handy. By the fact that NAd and N'M' are both pushouts, 
there exist unique morphisms a and r making the following diagrams commute: 




f M, (NM), 



ML 



{N'M'), 




(NM), 



We now use a standard trick of category theory to see that these unique a and r are isomor- 
phisms, and in fact inverses of each other. Using the uniqueness part of the pushout property 
we know the identity chain map (iVM),— ^(ATM), is the only morphism for which the diagram 




commutes. But the diagram also commutes if the dotted arrow is the map to a, and therefore 
Toa = 1. An identical argument shows that aor is the identity chain map {N' M'),-^{N' M'),. 
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We thus have an isomorphism a: (NM), (N'M'), such that the diagram 

S, ^ (NM), ^ 5^' 



(iV'M'). 

commutes. It is not hard hard to check that a preserves the inner product on p-cochains, so 
that we have NM ~ N'M' 

Now suppose we have three chain cobordisms 



M, 



s', s, m; s, — ^ Mi! s'^ 



Showing that composition is associative involves constructing an isomorphism 

((M"M')M)). ^ (M"(M'M)). 

such that the diagram 

((M"M')M)), 





(M"(M'M)), 

commutes, so that the two parenthesizations are equal as equivalence classes. In fact, a is 
just the obvious isomorphism. Checking the details is tedious but not difficult. In fact, the 
equivalence relation in the definition of chain cobordism was motivated entirely by the need 
for associativity. 

The identity morphism from S to S consists of two identity chain maps 

1 



S, 

which in expanded notation is the diagram 

< Sq ^ Si 



■s. 



s. 




* Sn 


-1 


1 




■ 




* Sn 






-1 


1 




< Sn 


-1 



Here we explicitly write in the position 5^ = in the middle row just to emphasize that the 
we are thinking of S both as an object, an (n — l)-complex, and as a morphism, an n-complex. 
We write this identity morphism as I5 : S — S. To see that this does satisfy the left and right 
identity axioms, suppose M: S ^ S', N: S ^ S". Then it is a routine exercise to check that 
the diagrams 

and 
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commute, where the vertical arrows are the obvious isomorphisms induced by the inclusions 
of S into M and N. 

This completes the proof that nChain is a category. A monoidal category is a category with 
a 'tensor product' so we must describe this product in nChain. In nCob the tensor product 
is just the disjoint union of spaces; in nChain the obvious analog is the direct sum of chain 
complexes, at both object and morphism levels. Given objects S and C, define S © C to be 
the chain complex S", © C,, with inner product on p-cochains given by 

(-,-) = (-,-)5 + (-,-)c. 

Similarly, for morphisms M: S ^ S' and N : C ^ C define 

M®N: S®C ^ S' ®C' 

to be the class of basis preserving and inner product preserving chain maps 

S,®C, — ^ M. © AT. ^ © C',. 

The identity object for © is the trivial chain complex 0,: 

0^ 0^ 

and the associator and unit laws are inherited from the obvious abelian group isomorphisms 
{G ® H) ® K ^ G ® {H ® K) SiuA G ® ^ = G ^ {) ® G. 

Making nChain a symmetric monoidal category requires that we also specify a 'symmetry' 
or 'braiding' which lets us switch the order in the tensor product. This comes from the abelian 
group isomorphisms G©-ff = -ff©Ginan obvious way. Checking that these constructions 
actually yield a symmetric monoidal category involves checking certain 'coherence laws' guar- 
anteeing the associator, unit laws and braiding get along appropriately. This a lengthy process, 
but not a difficult one. We refer the reader to [20] for the detailed definitions of monoidal and 
symmetric monoidal categories. ■ 

Theorem 2 Lattice p-form electromagnetism gives a chain field theory. 

Before proving this theorem, we remark that for cobordisms we must modify our definition 
of the action slightly. The idea of the modified action we shall use for cobordisms M : S ^ S' 
is that the action on p-cells in the 'boundary' S + S' should only count half as much as the 
action on cells in the interior. This keeps us from 'overcounting' the action when we compose 

two cobordisms S ^ S' ^ S" . liM: S ^ S' is a chain cobordism, let us rename our old 
naive action S, so that 

^^^^ = -^h^'n^i)n^3) 

in the M case. If A is a p-connection on the spacetime cobordism M: S ^ S', we then define 
the action on M to be 

S{A) = S{A) - ^-S{A\s) - ls{A\s,). 
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In the U{1) case, we don't define the action, but only its 'exponential' — the theta function 
in equation (4) in Section 7.3. 



g-S(A) ^ ^ ^-^h^^(A(dxi)-2nirc){A{dxj)-2njn) 

This leads us to define for a p-connection A on M: S S', 

g-S(A) ._ g-5(A)gi5(A|s)gi5(A|s,)_ 

where, ea just denotes the obvious thing — the reciprocal of the square root of the theta 
function "e"'^". 

Proof of Theorem 2: To define a functor 

Z : nChain^ Hilb 

we must specify what Z does to objects and morphisms. When space is the (n — l)-complex 
S, the classical configuration space of p-form electromagnetism with gauge group U{1) is the 
space of p-connections on S modulo gauge transformations: 

AjS) _ CP{S,Uil)) 

g{s) BP{s,u{^)) 

Quantization gives the Hilbert space of states, the space of all square-integrable functions on 
classical field configurations. We thus define 

Now Z must also assign to each morphism M: 5 — S" in nChain a linear map Z{M) : Z{S) — 
Z[S') between Hilbert spaces, corresponding to 'time evolution' of states on the slice S of space 
to the slice S' . To determine Z{M) it suffices to specify, for any ij; € Zi^S) = Lp'{A{S)/Q{S)) 
and 4> e Z{S') = L'^{A{S')/g{S')) the transition probability 

Ja(m) 

where for convenience we absorb the normalization factor into the definition of the measure 
VA — we can do this since the partition function converges when the gauge group is U{1). 
Also because we are using ^(1), factoring out gauge transformations is not necessary for the 
path integral to make sense. Hence we integrate over A{M) rather than A{M)/Q{M). Since 
U{1) is compact and the integrand is gauge invariant, the result is the same. 

To show functoriality we must show that this definition of Z respects identity morphisms 
and composition. First suppose we apply Z to the identity cobordism Ig - S ^ S, i.e. 

S, S, S, 
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Then given two states tpjcj) & Z{S) 



{A{S)/g{S)) we have 



Ja{s) 



I 4>{A\smA\s)VA = {<t>,i^) 
Ja{s) 



where we have used the fact that for the identity cobordism I5, we have 



g-5(A) ^ ^-S{A)JS{A)^lS{A) ^ ^ 



Hence Z(S') is the identity on Z{S) = {A{S) / G (S)) . 

As usual, showing that Z respects composition is the nontrivial part of showing Z is a 
functor. Suppose we have two composable morphisms in nChain: 



shows that a p-connection A° on the composite M'M is precisely a p-connection A on M 

together with a p-connection A' on M' such that A' and A agree when restricted to the 
common boundary S' . To enforce this agreement on S' in path integrals, we will use a delta 
function 6{A\s' - A'\s'). If {ijx : X £ A} is an orthonormal basis of L'^{A{S')/Q{S')), then 
expanding our delta function gives 



M: S' 
M':S'^ S". 



Since M'M is a pushout, the diagram 




6{A\ 



5' - A\s') = ^ {rjx{A'\s'),5{A\s> - A'\s')) m{A'\s') 



AeA 




^^a(^|s')'7a(^'|s') 



Now to show Z{M'M) = Z{M')Z{M), it suffices to show 



{(p,Z(M'M)ip) = {(p,Z(M')Z{M)ilj) 
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for all ^ G Z{S) and (j) G Z{S"). We have: 
{cj>,Z{M')Z{M)^l;) 

Ja{M') 

^ ^{A'y>)Y,{m,Z{M)^)m{A'\s')e-^^^'^VA' 
') AeA 

= {^,Z{M'M)^), 

as desired. Here we have used the expansion (7) and the fact that 

g-5(A)g-5(A') ^ g-S(A°)_ 

This is comes from the definition of the action on a cobordism in the case where the gauge 
group is M. Using this and the definition (6) of the inner product for a composite of cobordisms, 
we obtain: 

S{A) + S{A') = S{A) + S{A') - \ (s{A\s<) + ~S{A^\s')) \s{A\s) - \s{Al\s") 

= S{A°\m) + S{A°\m') - S{A°\s>) - ^S{A°\s) - 

= S{A°)-^S{A°\s)-\s{A°\s.) 
= S{A°). 

The result follows in the case where the gauge group is C/(l) from the relationship between 
Gaussians and theta functions. This completes the proof that Z is a functor. 

The essential features of a symmetric monoidal functor are a natural isomorphism 

$: Z{-)®Z{-)^Z{-®-) 

between the functors 

Z{-) (g) Z{-) : nChain x nChain ->■ Hilb 
and Z{— © — ) : nChain x nChain — Hilb 

and an isomorphism 

(f>- Iffilb ■^(OnChain) 

between the identity in Hilb and Z applied to the identity in nChain. The existence of is 
obvious: Z{0) is the space of functions on a one-point set, i.e. complex numbers, and C is 
the identity for the tensor product of Hilbert spaces. To complete the proof that the functor is 
symmetric monoidal, one must work out what $ is and check that $ and satisfy the required 
coherence laws. We leave these details to the reader, whom we again refer to Mac Lane's 
textbook [20] for the definitions. ■ 
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Appendix: Gaussian Integration 

The fundamental integral is the definite integral of the Gaussian function e~^^/^ over all of M: 



— oo 



More generally, if we consider the Gaussian e ^^/^^-^^ with standard deviation a, a simple 
change of variables shows that 



oo 



/oo 
-oo 



This gives us the usual normalizing constant for a one-dimensional Gaussian probability dis- 
tribution. Thus, using the properly normalized measure. 



27rcr 

we can calculate expected values of random variables. Given a function / of the Gaussian 
random variable x, the expected value of / is given by: 



y lira Jm. 



In general, such integrals are of course hard to carry out. Since any reasonably nice function 
can be approximated by polynomials, however, we will be content to calculate (,t"). It is not 
hard to show, using integration by parts and the fact that odd functions integrate to zero, that 

{x"") ^ ' ^n„-x'/2a' _ ] ^■''(ra - 1)!! n even 



\/27r(T J_oo I 



27rcr J-oo \ n odd 

Here, the double factorial (n — 1)!! means (ra — 1) • (n — 3) • • • (5) • (3) • (1), and we evidently 
need to define (—1)!! = 1 for the formula to be consistent with our original Gaussian integral. 

Next we generalize Gaussian integration to more than one dimension. The simplest case is 
the totally symmetric Gaussian bump, which has an obvious solution: 



/oo /"OO ^ 

. . . / e-2(^i+-+^n)dxi ...dxn = Y{\/2^= (27r) 
-oo J —oo 1 



fe=l 

which we can also write using the standard inner product on R" as 



Often, however, we want to do similar integrals where the argument of the exponential is 
some more general quadratic function of the .Xj. It is a basic result of linear algebra that any 
quadratic form q: F ^ M on an inner product space V may be written as 

q{x) = {x, Qx) 
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where Q is a self-adjoint linear operator V ^ V, i.e. a symmetric matrix^. The matrix Q is 
called positive whenever the quadratic form q is positive. 
We now calculate 



J 



where Q is a symmetric, positive definite n x n matrix. To do this, note that the positivity of 
Q implies Q has a unique positive square root — this is just the finite dimensional version of 
the square root lemma of functional analysis (of. Reed and Simon [26] p. 196). That is, there 
is a positive definite matrix B with B^ = Q. Moreover, since Q is symmetric, so is B, so we 
can write 

{x,Qx) = {x,BBx) = {B^x,Bx) = {Bx,Bx), 

where the adjoint B^ in the present case is simply the matrix transpose. To do the above 
integral, we can thus make the change of variables y = Bx, or in component form = BjX^. 
Using the proper Jacobian for the transformation, and the fact that det{B) = ^det{Q), we 
find 



jRn V det(Q) ■ 



To calculate expected values: If 

are m (not necessarily distinct) linear functionals on then we have the following simple 
formula for calculating the expected value of their product: 

{ili2...im) = \^^ [ hix)i2ix)...emix)e--2^''''^-^(tx 

= E n i^^Q'^'^j) 

pairings pairs £i,£j 

The proof of this relies on the facts that both sides of this expression are linear in each of the 
ii , and two multilinear functionals agree if and only if they agree when all of the linear factors 
are the same. 

In particular, for m = 2, if we take Q to be the identity on so that we are back to the 
symmetric Gaussian bump, we have a formula which, aside from the comma on the right side, 
looks like a tautology: 

{hh) = {h,h). 

For m = 4 we have 

(^lW4) = {£l,i2){i3,h) + {h,h){i2,h) + {h,h){i2,h) 



^In a more general setting, a self-adjoint linear operator is represented by a hermitian matrix, but since in this 
paper we are concerned with real vector spaces, the adjoint is merely the transpose. 
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where the 3 terms correspond to the 3 = (4 — 1)!! distinct pairings of £i, ... ,£4: 

Note that the number of distinct pairings of m objects is zero when m is odd and (m — 1)!! 
when m is even, so in the one-dimensional case where each of the linear functionals is just the 
identity map x and Q is just the number l/cr^, the formula reduces to 

/ m\ _ / — 1)!! n even 
I -\ n odd 

as we found before. 
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